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SUMMARY 


A package of computer programs has been developed which 
calculates the In-duct acoustic modes excited by a 
fan/stator stage operating at subsonic tip speed. The 
following three noise source mechanisms are Included: 

a) sound generated by the rotor blades Interacting 
with turbulence ingested into, or generated within, 
the Inlet duct; 

b) sound generated by the stator vanes Interacting 
with the turbulent wakes of the rotor blades; and 

c) sound generated by the stator vanes interacting 
with the velocity deficits in the mean wakes of the 
rotor blades. 

The first two interaction mechanisms generate random 
noise (although most of the resulting energy may, in some 
circumstances, be clustered around multiples of the blade 
passage frequency), while the third mechanism generates 
tonal noise at the blade passage frequency and its 
harmonics. 

The fan/stator stage is modeled as an ensemble of 
blades and vanes of zero camber and thickness enclosed 
within an infinite hard-walled annular duct. The acoustic 
pressure within the duct is calculated by distributing 
pressure dipoles on the surface of the rotor blades or 
stator vanes and calculating the pressure at an arbitrary 
point within the duct via the normal mode expansion of the 
Green's function for an annular duct. By this procedure one 
obtains an Infinite series for the sound pressure within the 
duct. Each term contains a normal mode of the duct 
multiplied by the amplitude of that mode. The amplitude (or 
the expected value of the amplitude spectral density) of 
each propagating mode is computed and summed to obtain the 
sound power flux (or the expected value of the sound power 
spectral density) within the duct. These calculations are 
carried through for both upstream and downstream propagating 
modes. 
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The equations relating the duct acoustic mode 
amplitudes to the pressure dipole distribution on the rotor 
or stator blades make no assumptions about the ratio of the 
wavelength of the sound generated to the blade chord. 

However, to simplify the computation of the dipole 
distribution generated by a given periodic or random inflow 
variation, the so-called strip theory approximation is 
used. Specifically, the dipole distribution is calculated 
by deleting, from the convected wave equation, all terms 
containing derivatives with respect to radius, but retaining 
the radius as a parameter in the boundary conditions (which 
derive from the blade geometry and the Incident fluid 
flow). At each radius, therefore, the equations to be 
solved are those for a linear cascade in subsonic flow. 
Although radial derivatives are deleted from the wave 
equation in calculating the chordwlse pressure distribution 
on the blades, radial variations in the amplitude and phase 
of the pressure distribution are taken into account in 
integrating over the blade/vane surface to obtain the 
amplitudes of the propagating duct modes. 

Turbulence drawn into or generated within the inlet 
duct is modeled as nonhomogeneous and anisotropic random 
fluid motion, superimposed upon a uniform axial mean flow, 
and convected with that flow. In the computer programs, the 
inlet turbulence velocity auto-correlation function is 
computed by a set of subprograms. It is Intended that the 
form of the correlation function be chosen to fit turbulence 
data collected from a transducer fixed in the inlet duct (as 
opposed to rotating with the fan). To allow some flexi- 
bility in selecting the functional form of the velocity auto- 
correlation function without requiring an excessive amount 
of data for its definition, the auto-correlation function 
has been expressed as the product of three functions. Each 
function depends upon only a single variable, one of the 
three cylindrical polar coordinates of a point within the 
duct: r, (Ji, x. Rational methods for selecting these three 
functions are suggested in a companion report (Ref. 1). 

The flow downstream of the rotor is also modeled as 
anisotropic and nonhomogeneous random fluid motion, 
superimposed upon a steady, but in this case spacially 
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nonunlfopm, mean flow. The mean (time-averaged) component 
of the flow Is assumed to have no radial component, so that 
the mean flow streamlines form cylindrical surfaces. It is 
further assumed that if one of these cylinders is unwrapped 
to form a plane, the mean flow streamlines will be 
parallel. This is equivalent to ignoring the effects of 
viscous and turbulent diffusion in the rotor wakes, a 
Justifiable approximation as regards the calculation of the 
forces on the stator vane, provided that the wake thickness 
chosen is that which is obtained at the axial station of 
the stator vanes themselves. The magnitude of the mean flow 
velocity is constant on lines parallel to the mean flow 
streamlines, but varies periodically in the direction normal 
to the mean flow streamlines. 

Rotor wake turbulence is modeled as anisotropic and 
inhomogeneous turbulence carried along with the mean flow 
downstream of the rotor. The intensity of the turbulence is 
constant on lines parallel to the mean flow streamlines, but 
periodically varies in the direction normal to them. The 
turbulence velocity auto-correlation function is expressed, 
in rotor-fixed coordinates, as the produc*’ of three 
functions, jach of which is a function of nly one of the 
three spatial variables r, 4> , x. Just as for inlet 
turbulence. 

Equations for the duct mode amplitudes, or expected 
values of the amplitudes, ar«-. derived in Volume 1, while 
documentation for the computer programs is supplied in 
Volume 2. 
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CHAPTER 1 
INTRODUCTION 

This is the final report of a project designed to 
develop a computer program which calculates the modal 
content of the acoustic fields set up in the inlet and 
exhaust ducts of a turbofan by the actions of nonsteady 
vortical velocity fluctuations convected past the rotor and 
stator blades. Such a program, supplemented by mathematical 
models of the impedance of the inlet and exhaust termina- 
tions, can be used to calculate the far field noise. The 
program can also be used to provide the initial conditions 
rea.uired for the investigation of duct liners, particularly 
liners tailored to suppress a specific mode or set of modes. 

The fluctuating flows come from three distinct 
sources: 1) turbulence drawn into the inlet duct, 

2) the turbulent component of the rotor blade wakes, 
and 3) the mean component of the rotor blade wakes. 

The first two motions are random, while the third is 
periodic (as seen by an observer fixed in the duct). 

It is characteristic of all three flows that the pressure 
fluctuations associated with them are an order of magni- 
ude smaller than the velocity fluctuations, and that all 
three are embedded in, or convected with, the mean flow. 

When the velocity fluctuations are carried past a rigid 
surface, such as a fan or stator blade, the requirement 
that the flow conform to the shape of the blade sets up 
additional (acoustic) disturbances, which propagate both 
upstream and downstream in the duct. 

The acoustic pressure in the duct can be calculated in 
a straightforward manner, using the Green's function for the 
duct, provided that the pressure distributions on the 
surfaces of the rotor or stator blades are known. If the 
Green' s function is expressed as an infinite series of the 
normal modes of the duct, the modal amplitudes of the 
acoustic pressure are obtained directly. The problem 
remaining, therefore, is to determine the pressure 
distributions on the surfaces of the rotor or stator blades 
generated by the convected vortical velocity disturbances. 
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This is accomplished by modeling the blades as surfaces of 
zero thickness and camber, which support an unknown but 
continuous distribution of pressure dipoles. The particle 
velocity normal to the blades on the blades themselves ■’s 
calculated from the unknown dipole distribution; this 
velocity is then required to nullify the normal comp'’".ent 
the vortical disturbance velocity (so that the total fluid 
velocity conforms to the shape of the blades). The result 
Is an integral equation for the dipole distribution, which, 
depending upon the approximations introduced, may admit e/n 
analytic solution, or may require numerical treatment. 

The approximations required to arrive at an analytic 
solution are considerable: essentially they ignore the 

annular or circular geometry of the fan (or stator), the 
three dimensionality of the disturbance, the interaction 
between blades, and the compressibility of the fluid. What 
remains is incompressible flow around an isolated airfoil 
in two dimensions, the solution of which, for a sinusoidal 
vortical gust, is known as Sear’s function. This was the 
approach used by Kemp and Sears (Refs. 2 and 3) in their 
original investigations of the interactions between blade 
rows in axial flow turbo-machinery. The incompressible 
approximation has since been shown to be unacceptable. 

For example. Fleeter (Ref. 4) found that compressibility 
can change the pressure on the rotor or stator values by 
as much as a factor of two, while Kaji (Ref. 5) calculated 
even greater changes (20 dB; a factor of 10) in the sound 
pressure level upstream of the blade row due to source non- 
compactness . 

At the opposite extreme, as regards complexity, are 
calculations carried out by Kobayashi (Ref. 6) and Koba- 
yashl and Groeneweg (Ref. 7), using equations derived by 
Namba (Refs. 8, 9) for an annular blade row in compres- 
sible flow. Because of the annular geometry the kernel 
function, which relates the dipole distribution to the 
velocity normal to the blades, is elaborate in form and 
time-consuming to compute. The associated integral 
equation is two-dimensional (i.e., both the spanwise and 
chordwise distributions of pressure dipoles must be 
determined). In his numerical work, Kobayashi made 
comparisons between his "exact'' aerodynamic theory and 
various approximate methods of calculating the normal 
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component of the Induced velocity. One of these, which he 
called the "quasl-three-dlmenslonal” approximation, 
coincides with what most would call a "strip theory" 
approximation. In calculating the Induced velocity, 
derivatives with respect to radius are deleted from the 
convected wave equation, so that one is left with the kernel 
function of a linear cascade In nonsteady compressible 
flow. The particular geometry of the cascade, as well as 
the vortical Inflow velocity, depends on the radius, and 
this parametric dependence is retained in calculating the 
Induced velocity (and hence ultimately the dipole 
distribution). The effects of radial or spanwls' variations 
In the amplitude and phase of the vortical disturbance 
velocity, as well as radial variations In blade chord, 
inter-blade gap, and stagger angle, are thereby taken into 
account, even if in an approximate manner. Kobayashl found 
that this procedure introduced errors no greater than 2 dB 
In the computed magnitudes of the acoustic modes set up In 
the duct upstream and downstream of the rotor. 

The strip theory approximation, as described above, has 
practical advantages. It is considerably easier than the 
"exact" method to implement numerically, and it accommodates 
more conveniently a variety of vertical inflow disturbances 
(e.g. random as well as deterministic motions). For these 
reasons, it has been selected for use in the computer codes 
developed in this project. 

The role of atmospheric turbulence in generating either 
broadband or tone-like noise, depending upon the extent to 
which the turbulence is distorted as it is drawn into the 
turbofan inlet, is by now widely known and well understood. 
Measurements made by Hanson (Ref. 10) and othe'^s have shown 
that when initially isotropic (or nearly isotropic) turbu- 
lence is drawn into a duct moving at zero or very low for- 
ward speed, the streamwlse component of the turbulence 
velocity is reduced as compared to the transverse compo- 
nents. Elongated turbulence eddies are intercepted by 
many rotor blades as they pass through the fan: because 
each blade experiences essentially the same velocity 
fluctuation, the noise generate"? is tone-like in nature, 
containing prominent peaks at multiples of the blade 
passage frequency. On the other hand, when the 
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forward speed of the turbofan is comparable to the velocity 
within the duct, the turbulence remains roughly Isotropic, 
and the spectrum of the sound generated Is less peaked. 

Modeling Inlet turbulence rp'^ulres that the velocity 
correlation function of the turbulence be anisotropic. We 
have chosen to write the turbulence velocity correlation 
function as the product of three functions, each of which 
depends only upon the quotient of one space coordinate and 
the corresponding length scale. This model Is not based 
upon or suggested by any physical model of fluid turbulence, 
but has the advantage of being easily adjusted to fit 
measured turbulence spectra, without requiring an excessive 
number of measurements for the purpose. Other models of 
inlet turbulence, such as the axlsymmetric turbulence model 
used by Kerschen (Ref. 11), could easily be Incorporated If 
desired. 

Laksmlnarayana et al. have collected an extensive set 
of measurements of the mean and fluctuating components of 
the flow downstream of a uiultlbladed fan operating at 
subsonic tip speed, using a transducer rotating at the same 
rate as the fan (Refs. 12, 13, 1^). Their measurements of 
the mean (or time-averaged) component of the flow at loca- 
tions remote from the blade hub or tip Indicate that the 
radial velocity immediately downstream of the fan blade 
trailing edges is substantial, but that this velocity 
decays rapidly, so that at distances greater than about 
one-half blade chord dovmstream the circumferential and 
axial components of the mean velocity predominate. Measure- 
ments of the velocity defect profiles show marked asymmetry 
about the streamlines on which the minimum velocity occurs, 
but this asymmetry also disappears within a half-chord length 
downstream of the fan blade trailing edges. At greater dis- 
tances, the normalized velocity defect profiles all show 
a Gaussian distribution, there being one such velocity de- 
fect profile for each fan blade. Due to the diffusion of 
momentum In the flow, the widths of the velocity defect 
profiles Increase v/lch distance downstream, while the 
maximum velocity defect decreases. But for the purpose -^f 
calculating the noise radiated by the stator. It Is suffi- 
cient to model the wake flow In the vicinity of the sta- 
tor. This can be done by ignoring momentum diffusion. 
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provided that the wake thlcl less and velocity defect are 
assl' 5 ned the values they attain in the vicinity of the 
stator, say at the leading edges of the stator vanes. 

With these approximations, the mean flow streamlines 
down-stream of the fan become a series of parallel lines, as 
viewed by a cylindrical surface opened out to form a flat 
plane. The velocity defect profiles are Gaussian, with the 
wake width and maximum velocity defect left as parameters to 
be specified. The flow velocity, being periodic in the 
azimuthal direction. Is conveniently expressed as a Fourier 
series for the purpose of calculating the inflow velocity 
fluctuations normal to the stator vanes. 

One additional aspect of the mean wake geometry is note- 
worthy. At each radius the mean flow streamlines all bear 
the same angle to the centerline of the duct, but this angle 
is a function of the radius, the precise nature of which 
depends upon the radial variation of design lift coefficient 
and chord length of the rotor blades. The locus of center- 
lines of any given wake form a surface whose shape depends 
upon the radial variation of the wake angle; this surface 
will in general intersect any selected stator vane at only 
one point, and this point will move along the leading edge 
of the stator vane as the fan rotates. Particularly if the 
gap between the fan and stator is large, each stator vane 
will at any given instant intersect many rotor blade wakes, 
and for each wake intersected, there will be two changes in 
the sign of the disturbance velocity seen by the stator 
vane. The pressure induced on the stator vane will have as 
many sign reversals as the inflow velocity, and will tend to 
excite duct acoustic modes with that same number of circular 
nodes. If, at the frequency of interest (which of course 
must be the blade passage frequency or one of its 
harmonics), few or none of the duct modes with that number 
of radial nodes propagate in the duct, the stator will not, 
at that frequency, be an effective generator of sound. The 
significance of this wake roll-up phenomenon in the genera- 
tion of sound by the stator was pointed out previously by 
Bliss, et al. (Ref. 15). In their discussion, they chose 
to emphasize the radial trace velocity of the vane/wake 
intercept instead of the number of nodes in the normal 
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component of the Inflow velocity. Their criteria for 
efficient sound generation, namely that the trace velocity 
of the points of Intersection of the wake centerlines with 
the stator vanes be supersonic, is perhaps more appropriate 
to an unshrouded rotor/stator combination, but the two 
criteria are roughly equivalent when the ratio of hub and 
tip radii is close to one, and are exactly equivalent for a 
set of two linear cascades in relative motion between two 
infinite parallel planes. 

In the equations derived in this report, and in the 
computer programs developed from them, the radial variation of the 
the angle between the mean wake streamlines and the duct 
centerlines may be specified arbitrarily, so that the 
effects of wake roll-up are properly accounted for. 

Wake turbulence measurements made by Laksmlnarayana, et al. 
(Rtfs. 13, 1^) using a transducer rotating at the same speed'” as the 
rotor indicated that the turbulence on the centerline of the 
mean velocity defect wakes is both more intense and contains 
more energy at high frequencies than does the turbulence in 
the region between wakes. The maximum turbulence Intensity 
coincides roughly with the maximum mean velocity defect, and 
the fall-off in intensity with distance normal to the wake 
centerline follows roughly the same Gaussian distribution as 
the mean wake velocity defect. The width of the turbulent 
wake is also roughly the same as the width of the mean wake. 

Based upon these observations, we have chosen to model 
rotor wake turbulence as homogeneous turbulence modulated by 
a function which is constant on the wake centerlines, but 
which varies periodically along lines drawn no1?mal to the 
wake centerline. This model is at variance with the 
observations quoted above in that the turbulence midway 
between wake centerlines has the same frequency content (but 
at diminished intensity) as the turbulence on the center- 
lines. But because of the diminished intensity of the 
"between wake" turbulence, it cannot figure significantly in 
the noise generated by the stator. 

Only one other model of rotor wake turbulence (as ; 
applied to acoustic calculations) has been published to - 
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date. In his model, Hanson (Ref. 16) took a very different 
tack by modeling both the mean and turbulent components Qf 
the wake flow as a series of velocity defects of randomly 
varying amplitude and position. The mean value of the 
amplitude of the wake defect coincides, of course, with that 
of the mean wake velocity defect, and the functional form 
(or shape) of the velocity defect must likewise conform to 
that of the mean wakes. In this regard, Hanson's model Is 
less flexible than the one used In this report. Also, In 
calculating the noise radiated by the stator, Hanson modeled 
the stator vanes as point dipoles. He was therefore unable 
to account for the (quite reasonable) possibility that the 
radial Integral scale of the turbulence might be consider- 
ably shorter than the span of the stator vanes. The point 
dipole model also requires that the wavelength of the 
radiated sound be considerably greater than the chord of the 
stator vanes, which at high frequencies Is not the case. 

No attempt has been made to model the flow field 
downstream of the rotor hub or tip; such work Is beyond 
the scope of this Investigation. 

The remainder of the report Is arranged as follows: 
pertinent aspects of the geometry of the duct, rotor and 
stator are alscussed In Chapter 2. In Chapter 3, a brief 
discussion of the acoustics of annular ducts Is presented, 
along with a derivation of the amplitudes (or expected 
values of the amplitudes) of the duct modes excited by a 
given pressure distribution on the fan or stator blades. 

The computation of the blade or vane pressure distribution 
caused by a given vortical inflow velocity Is outlined In 
Chapter and mathematical models of the fan blade mean 
velocity defect wakes, rotor inlet turbulence, and rotor 
wake turbulence are provided in Chapters 5, 6, and 7, 
respectively. Finally in Chapter 8, selected aspects of 
the numerical procedures used in the computer programs 
are discussed. Program listings and documentation are 
supplied in Volume 2. 
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CHAPTER 2 


GEOMETRY OF DUCT, ROTOR, AND STATOR 


For our purposes, a turbofan is a rotor and stator 
combination mounted in an annular duct of infinite length 
(i.e., reflection of acoustic waves from the ends of the 
duct are to be ignored). A sketch of the rotor/stator 
combination is shown in Pig. 1. As shown in the figure, a 
cylindrical polar coordinate system is established in the 
duct, with the polar axis lying along the duct- centerline. 
The axial coordinate x increases in the direction of air 
flow, and the rotor rotates in the sense of increasing 
The outer radius of the duct is rp and the inner radius 
is ry. The rotor has B identical evenly spaced blades, and 
the stator has V identical evenly spaced vanes. 

The rotor blades and stator vanes are modelled as 
twisted sheets of zero thickness and camber, whose pitch 
angle and chord vary with the radius, r. The axial and 
azimuthal sweep, if any, of the rotor blades are defined 
by extending a radial line from the axis of rotation 
through the mid-chord point of the blade at the hub 
(radius » . The location of the mid-chord point at 

any other radius r is defined by the two component vector 
g = (61,62) which gives the displacement of this point 
from the radial line. (61 = axial displacement and 
62 * azimuthal displacement.) In general, both 61 and 
62 are functions of the radius r, and by definition 
61 * 62 = 0 at the hub (r = Tt^) • See Figure 2. 

The intersections of the rotor blades and stator vanes 
with a cylindrical surface of radius r is shown in Fig. 3. 
The stagger angle of the blades is Xj and that of the vanes 
is 0 (both are functions of r) . The spacing between the 
blades is 2irr/B, and between the vanes, 2iTr/V. The local 
semichord of the blades is b (a function of r) . The same 
symbol is used for the semichord of the vanes. The dis- 
tance between the rotor and stator is d. 
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Figure 3. Rotor/Stator Geometry 
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CHAPTER 3 
DUCT ACOUSTICS 


3.1 Normal Modes in an Annular Duct 


The normal modes of an annular duct with hard walls are 
the set of solutions to the two-dimensional Helmholtz 
equation having the special form. 

yCrjif) « f(r) exp(ira((») (1) 


where m is any Integer, positive, negative, or zero. The 
Helmholtz equation is, in polar coordinates. 




+ 


i- 1L\ 

r^ 3(^2/ 


'i + = 0 


( 2 ) 


In this equation < is an undefined constant. As is shown 
below, the equation has non-trlvial solutions which fit 
the appropriate boundary conditions on the walls of the 
annular duct only for certain specific values of <, 

If Eq. (1) is substituted into Eq. (2), we obtain 


d^f + 1 ^ 
dr 2 r dr 



(3) 


Substituting u = Kr for r reduces Eq. (3) to Bessel’s equa- 
tion of order m: 



The solutions are mth order Bessel functions of the first and 
second kind, i.e., f(u) •» where 

- AJjjj(Kr) + BYj^(icr) . (4) 


If we specify that 'F(r,<j») • ixpClra*) is either the 

pressure or the velocity potential, then the radial deriva- 
tives of 4i-j must vanish at the inner and outer walls of the 
duct, so if A and B are to be not both zero, 
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= 0 


(5) 


This transcendental equation has a countably infinite number 
of roots K for every Integer m; if we denote these roots 
by <m.n» ^ “ 1,2,3,..., and arrange them in order of 
Incredsing magnitude, then the functions n^) have 

(n-1) zeros in the Interval rjj £ r £ rj). These functions 
are also orthogonal with respect to the weight function r 
over the same interval; that is, 

f^D 

■■h 


unless A a n. Using this fact, plus the fact that the func- 
tions exp(lm(|)), m = any integer, are orthogonal on the 
interval 0 £ 4» £ 2tt, it Is easy to show that the normal modes 
ipm(Km exp(lm4>) are orthogonal over the cross-section of 
the ddet. Specifically, 


J J ^ j^r)exp( -lk<J))}rd(|idr 

= 0 1fkj^morJ.?^n 



„r)r dr 
m,n 


( 6 ) 


if k » m and £ ■ n. It is convenient to adjust the constants 
A and B in Eq. (4) [which are determined by Eq. (5) only to 
within an arbitrary multiplicative factor] so that 


••h 



1 


The orthogonality relation [Eq. (6)] then becomes 


13 


ORIGINAL PAGE 13 
OF POOR QUALITY 









( 8 ) 


A special FORTRAN program has been written which cal- 
culates <ni n constants A and B In Ea. (^) subject 

to the norAallzatlon given In Eq. (8); for details see 
Appendix C. 


The significance of the normal modes is that they can 
be used to represent pressure patterns which propagate 
within the duct without change in form, and that any acoustic 
field within the duct, however generated, can be repre- 
sented as a suitable combination of these patterns. 

The wave equation in a duct containing a fluid moving 
at a uniform axial velocity U ls» 


7^p + 


ii£ 




p « 0 


where p is the acoustic pressure, and 7^ is the two- 
dimensional Laplace operator. 


( 9 ) 


5 ti-lL . 

Now assume that p is the real or imaginary part of 


V^'^m,n’^^ exp[l(m<|i-YX-ait)] 


( 10 ) 


This pressure pattern hasm diametral nodes (like the spokes 
of a wheel), and (n-1) concentric circular nodes. If we 
substitute expression (10) into Eq. (9), and recall that 

7*{i|;^(k _r)exp(im^)} - -<* „'*»„(< _r)exp(lm(|>) , 

m m,n m,n m m,n 


the following relation is obtained; 


•Between the rotor and stator an appreciable circumferential 
velocity exists as well. This velocity is not accounted for 
in Fig. ( 9 ). 
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Solving for y In terms of k and w, we find that 
Y = where 




( 11 ) 


Equation (11) shows that y is real whenever w is real and 
w/co > 3iCjjj yj. Under these conditions, the pressure pattern 
propagates ’unchanged along spiral paths normal to the lines 

mi}) ~ yx - lot = const . 

If to is real but lo/co < n» then y is complex and the 
pattern grows or decays exponentially along the duct depend- 
ing on which sign is selected in Eq. (11). 

Modal pressure patterns such as 

can be superimposed to form a general description of the 
acoustic pressure field in a turbomachine. Thus, if we 
multiply the pattern above by an arbitrary function of lo , 
say Pm n^^^)> integrate over frequency (io/2it), and sura 
over mjn, we obtain* 


* I ^ 

m n ’ ^ 


Pmn^“)®^P^it'^<J>‘‘'Tn jjj(w)x-u»t]} ^ (13) 


where now x is the point (r,x,i|)) in the duct. Note that in 
general p(x,t) is composed of both propagating and non- 
p’^opagatlng modal pressure patterns. The Fourier transform 
of this equation, namely. 


P(x,u) - I I ?na^(w)'(»j„(<jn ^r)expCim^+iy^ 
m n ’ * 


*Where no integration limits are specified, -» to +• is 
implied. 
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exhibits mere explicitly the role of the functions 
as complex modal amplitudes. 

3.2 Acoustic Energy Flux 

In a duct containing a fluid flowing at a uniform 
axial velocity U, the Instantaneous acoustic intensity 
is (Ref. 17 ) . 


I * 



+ Uu J( PpU+Up) 


(15) 


where Po is the nominal fluid density and p, p, and u 
are the Instantaneous acoustic perturbation pressure, 
density, and axial velocity, respectively. The acoustic 
pressure and density perturbations are proportional to 
one another (p = c^p), so Eq. (15) can be written in 
terms of the pressure and axial velocity only: 


I = (l+M^) pu + 


M 


Po°o 


+ p c Mu' 
*^0 0 


( 16 ) 


where M = U/c^ is the nominal axial flow Mach number. 

Now Introduce the Fourier transform s> of p and u, 
defined as follows: 


p(x,a>) 



p(x,t)exp(io)t )dt 


(17) 


p(x,t) *■ I p(x,u))exp(-lwt) (18) 

with the corresponding definition for u(x,a)), (In both 
integrals above the limits of integration are understood 
to be ±“. This convention will be adhered to throughout 
this report.) The axial intensity can now be written as 
a double inverse Fourier transform: 


I 


|J{ (1+M* )p(x,w)u*(x,v) 
+ PjCjM u(j 5 ,u)u*(x,v)> 


+ p(x,w)p(x,v) 

Po^O 

• expC-i(u)-v)t] 1^ 1“ . 


(19) 


i 


/ 

IT 
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where the superscrlp: ‘ denotes the complex conjugate. 
Both p(x,w) and'UCx,;'!; . m be written as sums of modal 
pressure patterns, ~a:i iSq. (14). Thus, 




U(x,v) .4^(<j^^j^r)exp[lk<^-iYj^^j^(v)x] 


The modal co-.-ff iclents Pmn^“) inde- 

pendent, but are related through the^axial momentum 
equation. 


3t ^ ^ 3X ^ p, 9x ^ * 


If Eqs. (20) and (21) are substituted into the Fourier trans- 
form of this equation, the following relationship is obtained: 


mn 


Tf P ( ^ ) 

pijU ^mn^ ^ 


where 


A (u) 5 


(23) 


The axial intensity can now be written entirely in terras 
Pim= 




(aH) 


d(u dv 


• expU[Y„_^(<o)-Vj_^(v)]A-l(<a-v)t]).p^(.,)pJj(v) 
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The Instantaneous sound power is found by integrating 
the intensity over the cross-section of the duct. That 
is: 


Power 



rd'i>dr 


7r(r^-r^ 


m n*' ^ 


CM^-(l+M^)A« (v) 
WJl 




( to ) -Y * 
n,m 'n,m 

dtp dv 
2 it 2 tt 


(v)]x 


( 25 ) 


The noise gen-*rated by the stator blades interacting 
with the mean wake momentum deficits of the rotor is 
periodic, the period being the blPde passage Interval. 
Thus 


p (to) = 2it y p 6(<o-nBs) 
^mn f; mns 

s 


= 2-.T I p»^^6(v-f2Br) 


mnr 


( 26 ) 


and the sound power flux in the duct is 

7T ( r ^ ^ ) 

Power « ^ ^ I J {M^-(l+M^ M* (rBJ^ ) 

^ m n r s 


+ A^(sKJ)AJ;_^(rB!!)]expti[Y^_^!5Bn)-Y»_JrB!l)]x 

-Ks-DBntllp^^pJ^p • (27) 


Those terns in this equation with r s represent fluc- 
tuations in the power. The time average sound power is 


Power ■ 


It ( r ^ — r ^ ) 




I I 

iu n 


l 


{m*-(i+mMa» (sBn)+|A„„(sBn)l ^)1 p_ J *. 

inn mn ' ‘mns* 

( 28 ) 
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The term in braces {•••} can be simplified: 

+M*6“(sBn/U)k „(sBn) 

{...}. JLiB (29) 


The upper set of signs apply upstream of the rotor (where 
the energy flows upstream or in the negative x-direction) 
and the lower set apply downstream. 

When the sound in the duct is random, Eq. (24) can be 
used to calculate the ensemble-averaged or expected value 
of the power, say <Power>, in terms of the_expected value 
of the mode spectrum in the duct, <Jmn(w) 

If the mode spectrum is correlated at differing frequencies 
(i.e., if <p;nn(<^) ^ ® when oj v), then the sound 

power within the duct fluctuates. For example, as will be 
shown later on, when turbulence is convected past either 
the rotor or the stator, a modal spectral density of the 
form 


<p (o))p* (v)> » 2Tr y p (co)5 (co-v-sBJi ) • 
^mn ^mn t mns 

s 


(30) 


is produced. If this expression is substituted into 
Eq. (24), the following result is obtained for the 
expected value of the acoustic power: 


TT(r2-r/J 


<Power> » \ I I exp(lsBnt)* [M*-(l-M* )A^(co+sBSl ) 

P A ^ m n Q 


mns 


(31) 

All terms with s 0 correspond to fluctuations in the 
power in the duct. If we average over one period, only 
the s » 0 term remains: 


<Power> 






do) 
?? • 


( 32 ) 


r 
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The integrand is the p^wer spectral density. That Is, 


:| 


I 

- t 


<f*ower> 


S (w) 


dec 

ri 


where 


S(u)) 






m n 


mno 


(w) . 
(33) 


But Pmno^***) Just the time average of the amplitude 
spectral density of the (m,n)th duct mode, which later on 
we will denote by Pmn(“)* Thus, we can write 


IT ( r* ^ ^ ) 

SM • I I {M*-(1-M2)A* (u) 

^0 m n 




(31) 


Again, the quantity in braces in this equation can be 
simplified: 


{ } 


+M26'*(w/U)k ^(uj) 


(35) 


3.3 Duct Acoustic Nodes Excited by Fluctuating Loads 

on the Rotor or Stator 

In this section, we shall derive equations which re- 
late the acoustic modes in the inlet or exhaust ducts to 
fluctuating loads on the rotor or stator blades. The 
starting point is Eq. (4.13) in Goldstein’s book A«ro- 
aaoustioa* (Ref. 17). 

*(3oldstein*s equation has the density fluctuation p(jc,t) 
on the left-hand side; the version shown on the following 
page is actually his equation multiplied through by cj. 


> 

; / 
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p(x,t) » [ f 7G(x,y,t-x) *f (y,x)ds(y) dx , ( 36 ) 

“m j J mm “W 

-« S ( X ) 


which gives the pressure fluctuation p(x,t) at any point x 
within the duct in terms of the force/unit area C exerted 
by the blades on the fluid. In this equation, G(x,y,t-x) 
is the time-dependent Green's function for a hard-walled 
annular duct. If we Introduce a cylindrical coordinate 
system such that the field point x is (r,Xi,(^) while the 
integration point y is (r' ,yi,$),~then G can be written 
as follows: 


G(x,y,t-x) 


571 ^ ^ 
m n 


„r)i|j (ic r ' ) 

m m^n m 


m^n 


exp[lm((|)-<|)) ] 




•f exp[-iw(t-x)+Yn , (37) 

J ^n,m 

where 


'n,m 


= 

' 8'c, 


k 

6^ 


sgn (y^-x, ) 


(38) 



(39) 


In Eq. (36), n^^ exp(im(|)) is a normal mode of the 

duct, and ipm(<m n^^^ exp(-lm<t»') is the complex conjugate 
of the same norAal mode, evaluated at (r',$ ) instead 
of (r,<|)). 

If the fluid is assumed to be inviscld, the force f 
will be normal to the surface. 


f - pn , 

t 
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where p is the local pressure. Thus, 

p(x,t) « II n(y) •VGCxjy ,t-x) p(y,T)dS(y) dx . (lio) 

S(y) 

The integration is to be carried out over both faces of 
each blade. On each blade, we can divide the integral 
in the equation above into two parts, one over the forward 
or upstream face of the blade, and the other over the 
downstream face. Calling the upstream face + and the 
downstream face we have, because the blades are very 
thin. 


P+n+ + P_l}i * (P+-P_) !} > 

where n is a unit normal vector erected on the middle 
surface of the blade. Define Ap as Ap = P_-P^. Then 

P(j»t) ““II n(x)*7G(x,x,t-x)‘Ap(x,x)dS(x)dx , (^1) 


where now the integration is over the blade surfaces (as 
opposed to the two blade faces). 

Because we are interested in frequency spectra rather 
than pressure time histories, it is convenient to calculate 
the Fourier transform of Eq. (4l), using the Fourier trans- 
form pair as defined in Eqs. (17 18). Before this can 
be done, however, we must first introduce the following 
chancre of variable: <|>' = $-ftx. The angle remains 

fixed as the rotor rotates, so the order of the time and 
space integrations in Eq. (kl) can be interchanged. 

With the above substitution, G contains x in the combina- 
tion expCi(u-mn)x] , so, carrying out the x integration, 
we identify 


Ap(;^* »T )exp[l(a)-mSl )x ]dx * Ap(^*,a'-mn) . (42) 
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Because t appears in Eq. ( 37 ) only as the a actor 
exp(-lwt), the (u integration is effectively an Inverse 
Fourier transform. Thus, 




H3) 


1 f 
2irk 

n.m 




]Ap(y,w-ntf2)dS(y) 


Note that 

„r)exp(im<j)-iY„ ) 
m m,n n,m i 


is a rotating pressure pattern of the type discussed pre- 
viously.* Equation (43) is thus a normal mode expansion 
of the Fourier transform of the acoustic pressure within 
the duct : 


Pmn^“^ ~ 2irk 


n,ra ' 


_r’ ) 
m ni>n 


m 


n(y)*7[exp(-im<j)'+iY^ jjjyi ) ]Ap(y,(o-mfl )ds(y) , (44) 


and 


p(;,u) - I I f)exp(lm4.-lY X,) . ( 45 ) 

m n ^ * 


*£xcept for the factor exp(-lwt). 
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In Eq. (^^), the integration is to be carried out 
over all B rotor blades. To reduce the region of integra- 
tion to one blade, arbitrarily select one blade as a 
reference blade. Assign this blade the number 0, and the 
remaining blades, in the direction of decreasing (fi*, the 
numbers 1 through (B-l)_j_ Let the pressure on the sth 
blade, s * 1, B-1, be APs(yo »w-mJ2) , where yo is 
the. point (r * ,y 1 , (j) • ) . Then the corresponding point on blade 
s is (r* ,yi ,(^’-2 tts/B. The total contribution of all B blades 
is thus 


Pmn^“^ “ 2irk, 




)-''Cexp(-W+lY„_^y, )]• 


fl^P.Cy. >4i-m(!)exp(l2iiins/B))dS(y ) . 

s s ..0 


(46) 


The normal vector on the reference rotor blade is 
n » (0,-sinx,-coox) , so 


n‘V[exp(-im(|)’+iYj^^jjjyj ] 


cosx-iYj^^jjjSinx)exp(-im(})’+iYj^^jjjyj ) 


To facilitate the integration over the blade surface, 
establish an intrinsic chordwise coordinate z’, varying 
from z’ « -b at the leading edge to z » +b at the trailing 
edge. Then, following the discussion of blade geometry 
in Chapter 1, coordinates yi and (j>’ are given by 

y^ = (S{ + z’cosx' 

(^* > -(6 *+z ’ sinx)/r • , (48) 

where x'(r') is the angle between the local blade chord 
and the axis of rotation, and (6^,6p is the displacement 









i 

I 
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of the blade mid-chord line from a radial line 
drawn through the mid-chord of the blade at the hub. 
Using r* and z* as Integration variables, Sq. (46) 
becomes 


d 


n,m ^ * 

n 


-b 


exp[i(Y^^j^cosx’ + ^ slnx ’ )z ' ]dz *dr 


(49) 


The corresponding equation for the sound generated 
by the stator is obtained from Eq. (49) by changing B to V, 
and by setting JJ = 0 and x' " 


■ 2fif 




n,m i 


m 


(<_ ,,r’ ) (— COS0 * 

m,ii \r 


■H 


+Y„ _slne ’ )exp(lY_ „5|+ —r 5 I)* 
'n,m *n,m i r ’ 2 


»+b 


-b 


V-1 

I Ap (r',z',o))exp 
s-0 ® 


12Trms 

V 


expCl(Y„ „cose’ - ^ sine ' )z' ]dz'dr' 
n r ’ 


(50) 
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ROTOR/STATOR CASCADE RESPONSE 


The pressure distribution on the rotor or stator 
blades can be calculated by dividing the rotor or stator 
into a series of radial "slices,” and calculating the 
pressure on each "slice" as though it were a linear cas- 
cade of thin flat plates. In this approximation, the 
radial variation of the Inflow to the blades is Ignored - 
an application of "sti io-theory" to rotor aerodynamics. 


The fluid velocity relative to the rotor blades (or 
stator vanes) is supposed to be a random or deterministic 
fluctuation u(r,xi,X 2 ,t) superimposed on a nominal flow 
Ur(r), which'is a function of the radius only. We are 
interested only in the fluctuating portion of the loading 
on the blades, which in the linear approximation is 
dependent only on the component of the velocity fluctua- 
tion that is normal to the blades, say w(r,x,t), where 
X = (xi,X 2 ). This velocity component can be expressed as 
a sum of traveling waves; 


w(r,x,t ) 


f r F 


w (r,k,w)exp[l(lc*x-u)t)] 


d^k do) 
(2ir)^ 2Tr 


(51) 


2 

where d k = dkidk2 

In the coordinate system shown in Fig. 4, the mid- 
chord points of the rotor blades are located at the points 
(0,sh), where s is any integer. The normal velocity that 
would exist if the blades were not present is found by 
substituting x « zc + sh into Eq. (51) > where c = (cosx>sinx) 
is a unit vector directed along the blade chords, and 
h = (0,h) is the vector separation between any two neigh- 
boring blades. The coordinate z is equal to -b at the 
blade leading edges and +b at the trailing edges. Thus 
the velocity that would exist if the blades were not pre- 
sent is 


w(r ,zc+sh,t ) 


exp(lsk*h) 


( fw(r,k,co)exp[i(k»cz-(ot )] , 

IJ " " *• (2tt)2 2^ 

(52) 


Because the blades are Impermeable, they Induce an 
additional velocity which negates the normal inflow velocity 
at the blades themselves. This induced velocity field has 
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associated with it a pressure distribution which is con- 
tinuous except on the blades where a discontinuity, say 
Aps(r,z,t), exists. The problem cf calculating the chord- 
wise pressure discontinuity (or pressure loading) on a 
cascade of blades exposed to the normal Inflow velocity 


w(r,k,co) expCik* (cz+sh)-lwt ] 

can be reduced to the solution of, an Integral equation 
of the following form: 


PjUj, w(r,k,o))exp[ik* (cz+sh) ] 


■■I’b _ 

K^(z-y) b 

-b 


(53) 


The kernel function K^, is derived in Appendix B. 

For our purposes, it is convenient to normalize Eq. (53) 

/V 

by dividing both sides by p U,„w(r,k,a))exp(isk*h) : 

0 r ^ 


exp(ik«cz) 


f+b (Ap^(r,y) exp(-lsk-h) ) ^ 

K^(z-y) I 2 

-b ( PpU^w(r,k,u)) ) 


(5^*) 


By assumption (see Appendix B) the quantity in braces is 
independent of s; we denote it by f (r,z,k,w). The solution 
of Eq. (5^) for f(r,z,k,u)) is obtained numerically. The 
chordwise pressure distribution is then 

— ' £. 

APg(r,z,k,w) * PjU^ w(r,k,(o)f (r,z,k,oj)exp(isk«h) . ( 55 ) 


Because the problem is assumed to be linear, the pressure 
loading due to an arbitrary normal Inflow velocity [as 
in Eq. ( 51 )] can be obtained by multiplying Eq. (55) by 
exp(-iwt) and integrating over w, ki, and ka. 

Thus, 


APg(r,z,t ) 



r,k,a))f(r,z,k,u)*exp(-ia)t+lsk«h) 
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/V 

/N 

Given the complex wavenumber-frequency spectrum w(r,](S, w), 

Eq. (56) allows us to calculate the pressure loading on each 

blade as a function of position (r,z) and time t. Where the 

Fourier transform (in time) of Ap is required, it is given by 

s 


\ 


^p. 


(r,z,o)) » Po^r If '^(r,k,(i))f (r,z,k,o)) •exp(isk»h 


h) 


d^k 
(27T) = 


(§7) 


For future reference, it is worthwhile to note here that 

w(r,k,u) is the triple Fourier transform of the Inflow to 
the rotor or stator, calculated in a coordinate system 
fixed to the rotor or stator, as the case may be. That is. 


A 


w(r,k,w) 


w(r,x,t )exp(la)t-ik* x)d^xdt 


(58) 


Because the inflow is necessarily periodic in x^, 
w(r,x,t) can also be written as a Fourier series; 


r /nx^v 

w(r,x,t) = I w^(r,Xj ,t)exp 


(59) 


where 


1 r 

w„(r,x, ,t) = w(r,x,t)exp j 

-irr 


dx . ( 60 ) 


If Eq. ( 59 ) is substituted into Eq. (58), we obtain 
£ ^ 

w(r,k,w) ■ 2 u J w^(r,kj,oj) 6 (kj - ~) 
and, in turn. 


n 
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Ap,(r,z,a)) 

o 


» p u 

*^o r 


^exp (isnh/r ) 
n 


w^(r,kj,(D) 


dkj 

f[r,z,kj,(n/r),w] ^ 


( 61 ) 


This expression for the pressure on blade s is convenient 
for calculating the duct modes excited by the rotor mean 
velocity deficit wakes, because it displays clearly the 
periodic nature of the rotor wake. 


i 
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MEAN ROTOR WAKE 


Of foo^. Q- 


, r 


5.1 Rotor/Wake Model 

In this section, we derive an expression for the 
complex amplitudes of the duct raodes excited by the Inter- 
action of the rotor blade wakes with the stator vanes. 

We are concerned here only with the mean value (time 
average) of the rotor wakes, which, being periodic, 
generate sound at harmonics of the blade passage rate, RB. 
Broadband noise generated by the turbulence contained 
within the rotor blade wakes Is dealt with In Chapter 6. 

Consider an Imagl ^ary cylinder of radius r, centered 
on the duct axis of symmetry. The Intersection of this 
surface with the rotor and stator Is depicted In Fig. 5. 

In this figure, two sets of coordinate axes are shown; 
axes (XijXa) are fixed In the rotor, while axes (xi,X 2 ) 
are fixed In the stator. The relation between them Is 


X » X + D + Jirte, 


( 62 ) 


where e.ls a unit vector In the azimuthal direction and 
D = D(r) Is the vector distance from a rotor blade to a 
stator vane, both selected arbitrarily. P Is defined 
only to within an additive vector r<l)e 2 , where repre- 
sents an arbitrary angle of rotatlon^of the rotor. 

So as to arrive at a reasonably simple model of the 
rotor wake, which will be valid In the vicinity of the 
stator vanes. It Is convenient to Introduce two plausible 
assumptions. They are 


1. no radial flow occurs, and 

2. pressure gradients and turbulent or viscous diffu- 
sion can be neglected over the chord of the stator. 


The fluid velocity then has only two components, 
(0,Wj,Wji In (r,X,,X 2 ) coordinates, and the equations of 
motion of the fluid reduce to 


(W‘7)W » 0 (63) 

7 • (p,W) • 0 . (6U) 


31 




A parallel flow of the form 




W = W(X*n)w , (65) 

v^here w and n are constant orthogonal unit vectors (see 
Fig. 5^) is a solution to these equations, and is the form 
we will assume for the wake. The velocity is everywhere 
parallel to the unit vector w = (cosx,sinx), and is con- 
stant in magnitude on the lines X*n = constant. The 
variation of the magnitude across the waKes is determined 
by the as-yet-unspecified function VI(X»n). 

W(X*n) is periodic, because all the wakes are assumed 
to be Identical, the period being the normal distance 
between wake centerlines, hcosx* W(X»n) can therefore 
be written as a Fourier series, i.e.^ 

W(X-n) - I ?■?) ’ 

where we have substituted 2irr/B for h, the gap between 
rotor blades, measured azimuthally. The unit vector n 
has components (-sinx,cosx) » so the Fourier series may 
also be cast in the following form. 


W 


I Wq exp 

q ^ 


(X^-X^tanx) 


(67). 


which indicates that the rotor wake velocity is periodic 
in the azimuthal airection also (as opposed to the direc- 
tion normal to the wakes). The azimuthal period is, of 
course, 2irr/B. 


In principle, the Fourier coefficients Wq of the 
wake flow could be determined experimentally by process- 
ing sufficient number of flow measurements collected 
from a transducer mounted behind the rotor, and if only 
the first few harmonics of the blade passage frequency 
are of Interest (it will transpire that the qth harmonic 
in the wake flow interacts with the stator to gene- 
rate sound in the qth harmonic of the blade passage 
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I'reAuenoy , this 1$ a feasible proposition. However, 
by assuminc; an e.xplicit land plausible^ form for the 
wake veioolty profile, it is possible to oharaoterise 
the wake velooiiy profile by only two independent at- 
tributes of the flow, as, for example, the mean wake 
velooily and the root-mean-square deviation from the 
mean velocity, or the mean velocity and the velocity 
deficit on the centerlines of the wakes. To do this, 
we express the wake velocity deficit as a series of 
"haystack" functions, all of the same form but dis- 
placed from onu anoth**!' by the gap between the rotor 
blades, h - Jnr/B. That is, momentarily setting X\ » 0, 
we write 


w(Xjl - A[1 - cl f lX,.-2TTrm,'Bl ] , (68) 

m 

where the constants A and c are yet to be evaluated, and 
the function f(X,) is a function such as 


fvX^^ ■ exp 6)*inJ] 




or any siml'.ar function whicli has the value 1 wheti 
■ 0, and wl'-lch decreases mottotonlcally to 0 as 
X) •* In Evj. lO'J), >5 is the ha If-velocity width of 
the mean wake, that is 


ft 


1 : 


t:o) 


If Eq. (o3) above is written as a Fourier series, 
we find that the Fourier coef I'lc lents are proportional 
to the Fourier transform of fCX,"), i.e,, 



q I* 0 


(71) 


The aeroth-order Fourier coefficient need not be cal- 
culated because, as with any Fourier series, it is equal 
to the mean value of the wake velocity, W. The constant 
eA can be evaluated either of two ways; one is to 
compute the root-mean-square value of the deviation of the 
wake velocity from its mean value. 
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Denote by Sj the sum 



and by (AW)j.nis the square root of the average of the square 
of the deviation of the wake velocity from its mean value. 
Then, the Fourier' coefficients of the wake velocity have 
the following values: 


W - TJ 


W. 


(AW) -s „ 

— ^ f(|S) , q ^ 0 


(72) 


and the velocity in the wake at any point (Xj,Xj) is 


W(X) - 

1 q^O 


?(|S) 


exp 


^ (Xj-Xjtanx) 


(73) 


An alternative scheme is to define eA in terms of the mean 
velocity W and the velocity at the wake centerlines, say 
Wq, Then, we find 


W(X) 


(W -W) ^ Q 

W + -4 T f(S^' exp 




qj»o 


^ (Xj-Xjtanx) 


(74) 
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where 



( 75 ) 


For the Gaussian profile of Eq. ( 69 ), whose Fourier trans- 
form is 


f(a) » 6 /Tr/ln2 exp[-(a6) /4 ln2] (76) 

the sums S, and S, are defined as follows: 

1 2 


Sj = 26/i7TiT7 I exp[-(mB6/r' V 4 ln2] (77) 

m*l 


I 2tT r 2 

S 2 J^exp[-toB6/r) /2 ln2]J . (78) 


5.2 Duct Modes Generated by Mean Wake/Stator Interaction 

To calculate the duct modes excited by the inter- 
action of the stator vanes with the mean rotor wake, we 
have first to calculate the loading on the stator vanes 
themselves. To do this, calculate the wake velocity 
relative to the stator vanes. 

Up = W - Grej , (79) 

find the component of this velocity that is normal to 
the vanes, 


w « U ‘N , (80) 

and substitute into this expression the coordinate trans- 
formation given in Eq. (62). To calculate the normal 
velocity on any specific vane, say vane s, substitute for 
X, 

•10 
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X * zc + sH . (8l) 

where again z is a chordwlse variable, equal to -b at the 
vane leading edges, and +b at the trailing edges. The 
result of these several steps is 

w = I Wq exp|ikqZ + - iqB^t j , (82) 

where 

Wq * Wq sin(e+x) exp 

and 

kq = ^ (sln0+cos0tanx) . (8^*) 


“ ^ (Dj-Djtanx) 


( 83 ) 


In these equations, D and are both functions of the 
radius r. 

The vane pressure loading generated by this normal 
Inflow velocity is obtained numerically by solving an 
Integral equation Just as discussed in Chap. 4. Speci- 
fically, let fq(r,z) be the solution of the following 
equat ion : 


exp(ikqZ) 



( 85 ) 


where Kg is the cascade kernel function derived in 
Appendix B. The parameters needed to calculate Kc for 
this situation are listed below: 

Frequency: w * qBK 

Chordwlse wavenumber: (sln0+tanxcos0) 

Inter-blade phase angle: a • 27rqB/V (See Appendix B) 
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Inter-blade gap: h « 2irr/V 

Stagger angle: -0 

Given the elemental stator vane loading function 
fq(r,z), the pressure loading on stator vane s is 


Apg ■ pjUj, Wq fq(r,z) exp(i2TtsqB/V) . (86) 


This is the qth harmonic of the pressure on vane s. To 
calculate the (m,n)th duct mode generated by this vane 
pressure, multiply both sides of Eq. (86) by exp(-iqBGt), 
calculate the Fourier transform, substitute this into 
Eq, (50), and invert back into the time domain. The 
sum over vane number in Eq, (50) can be evaluated 
explicitly; 


I exp|i^ (m+qB)| 

« V if ro + qB « pV 
» 0 otherwise, 

where p is any Integer, The final result is 

PoV Pd 

2%,n,q I ^ <3 

(5 cose+Y- ^ -Sine) expdv- _ +S62) (87) 

r ’n,m,q n,ni,q i v ^ 

r+b 

• J fq(i‘,2) ®xp[i(Yj^^^^qCose-“ sin0)z] dzdr , 

-b 

where 
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^n,m,q " ^j 
( = 1- (mm ± ic \ 

^n»m,q ' g2 \ c„ * *^n,ra,qj 


( 88 ) 

(8S) 


The quantity Pm n q complex amplitude of the 

(m,n)th mode at*tAe qth harmonic of the blade passage 
frequency. As is evident from Eq. (87) above. Pm,n,q is 
proportional to Wq, the qth-order Fourier coefficient of 
the mean rotor wake. In its present form, the computer 
code computes these coefficients from Eq. (74), which 
assumes a Gaussian profile for the individual rotor 
wakes, and which relates each harmonic in the Fourier 
series to the mean wake velocity and the velocity at the 
centerlines of the wakes. As mentioned previously, how- 
ever, the code can easily be modified to accept the Fourier 
coefficients of the wake velocity as inputs, and compute 
the duct modes directly from them. 


5.3 Sound Power 


The flux of sound power in 
tuting Eq, ( 87 ) into Eq. (28). 


the duct is found by substi- 
The result is given below: 


Power ■ 




z 

p 


E E 
n q 




mnq 


2 


where G(qBG) is the quantity in braces in Eq. (29), but with 
the integer s replaced by q. the equation above, the 
index m is related to the summation indices p and q through 
the equation m"pV-qB. Finally, if the sound power flux at 
one specific harmonic of the blade passage frequency is 
desired, the summation over q is deleted, and q is set equal 
to the desired harmonic number. 
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6.1 Inlet Turbulence Velocity Correlation Function 

Figure 6 shows a cyiindrical slice of the rotor which 
has been opened out to form a linear cascade. Axes (Xi,X2) 
are fixed in the duct, while axes (Xi,X2) move with the 
rotor blades. The equation relating the two sets of axes 
is 


X = X + nrte, . (90) 

The nominal Inflow velocity is assumed to be purely axial, 
so it can be written as U = Uei, where ei is a unit vector 
along the Xx and Xi axesT 

Let w be the component of the inflow turbulence which 
is normal to the rotor blades. We assume that this tur- 
bulence is convected with the mean flow y, l.e., a sample 
velocity field expressed in duct-fixed coordinates would 
be 


w(r,x-Ut) 

Many such samples are recorded and averaged to form the 
inflow velocity correlation function: 


<w(r^,x-yt) w(r2,y-Ut)> 




( 91 ) 


If we assume that the inflow turbulence is homogeneous at 
any radius r, and stationary in time, then 4 >d does not 
depend on x-Ut '.nd y-Ux separately, but only on their 
difference'*x-y- J(t-T)f we can write <!' as follows: 


*D “ ♦pCx-y-UCt-x) ,r,Ar ) 


(92) 


where 

r 

Ar 


|(r,+r^) 


r - r 
1 2 


ilO 


( 93 ) 
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The mean-square turbulence intensity is obtained by setting 
X = y and t = t. If is normalized so that 
then C]) is the root mean-square turbulence Intensity ex- 
pressed as a fraction of the nominal flow velocity in the 
duct: 


Aw * > 
‘ -TT- 


(94) 


6.2 Wavenumber- Frequency Spectrum in Rotor-Fixed Coordinates 


In order to calculate the power spectral density of the 
sound generated by the interaction of inlet turbulence with 
the moving rotor blades, we must calculate the expected 
value of the inflow wavenumber-frequency spectrum in rotor- 
fixed coordinates: 

A ft 

<w(r,,k,o)) w*(r.K,v)> = 

X Z ^ 

WWW W(r^,y-Ux)> 


exp(iujt-lvT-ik*X+iK*p d^X d^Y dxdt . (95) 

The expected value of the inlet turbulence is given in 
Eq. (9l)* li* we substitute this expression into the equation 
given above and Introduce new integration variables x and y, 
where 


X = X + nrte, 

Y = y + fine , (96' 


following Integral: 
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<w(r^,k,o)) w(r^,K,v)> » 

III III 

exp{l (u-Qrkj ) t - Kv-flrK^)! 

- ik*x + lK*y} d^x d*y dxdt . (97) 

Now substitute C * J 5 - y • y(t-T) for x, and integrate with 
respect to t, t 7 ano "recalling the following integral 
expression for the delta function: 

[ exp(lax) dx » 2 tt 6(a) . (98) 


The result is, 

1 ^ 

<w(r^,k,u) w(r2»K,v)> = 

(2TT)“eJU2 6(u)-firk2-y*k) 

/V 

• aCv-arKj-U-k) 6(K-k) Jp(k,r,Ar) , (99) 


where *D(lS»r,Ar) is the double Fourier transform of the tur- 
bulence velocity correlation function: 

Jj^(k,r,Ar) s j| ♦p(x,r,Ar) • 

exp(-ik*y) d*y . (100) 


Because the inlet duct is annular, the turbulence 
velocity correlation function is a periodic function, the 
period being 2irrfi2. can be written as a Fourier Series , 
l.e. , ^ 

♦jj(x,r,Ar) - ♦jjj(xi,r,Ar) exp (Imxa/r) ( 101 ) 


whei i 


4„(xi,r.4r) Sjij 


+Trr 


♦jj(x,r,Ar) 


-nr 

exp(-lmx2/r )dxj . 
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Now let the nonperiodic function <Kx,r,Ar) be defined as 
follows: 

(|)(x,r,Ar) = ♦jjCxjPjAr) if |x2l<Trr, (103) 

* 0 otherwise 

Then the Fourier coefficients of [Eq. (102)] can be written 
as the Fourier transform of ()>: 

(|)jj^(xi ,r,Ar; » ^$(xi ,m/r ,r ,Ar) (104) 

where ^ r 

(|)(xx,m/r,r,Ar) * (|)(x,r,Ar) exp(-irax2/r) dX2. (105) 

J ^ 

The double Fourier transform of is then 

r) • I I i(kj,s/r,r,Ar) 6(kj-s/r) 

** s 


Equation (99), the expected value of the complex 
wavenumber- frequency spectrum, contains five delta functions. 
[Note that 6(k-K) is a short-hand notation for the product 
of two delta functions: 6(k,-Kj) 6(k2-K2).l The five 

delta functions in Eq. ( 99 ) can be rewritten in the fol- 
lowing exactly equivalent form: 


6 (w-sfl-k'U) 6(v-w) * 


• 6(k-K) 6(k2-s/r) . 

The final expression for the complex 
spectrum is as follows: 


<w(rj,k,w) w*(rj,K,v)> 

€*U^2ir)'* 

d(v-u)) 


wav e numb er- f r e que nc y 


6(k-K) I 6(kj-s/r) ♦(k^ ,s/r,r ,Ar) 6(u)'Sn-U'k) (107) 

” s 


6.3 Inflow Velocity Spectral Density In Fixed And Rotating 
Reference Frames* 

Equation (107), which has been cast in the form most 
convenient for calculating the power spectrum of the sound 


•This section is adapted from Ref. 1, pages 37-47. 
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radiated by the rotating rotor blades, does not itself have 
an obvious physical Interpretation. It is worthwhile, 
therefore, to calculate the spectral density of the Inflow 
turbulence as "seen” by both fixed and rotating sensors. 
Comparing the two spectra provides some insight into the 
effects of rotation. 

The mean-square velocity of the inlet turbulence can 
be calculated by setting n » ra = r in Eq. (107), multi- 
plying it by (2ir)”® exp[i(v-u))t+i(k-K) «x] and integrating 
over all values of v, k, and K: 

^ m,0 * 


<w^(x,t)> = 

= Hill <w(r,k,w) w»(r,K,v) 


> • 


( 108 ) 


exp[i(v-uj)t+i(k-K) *x] - - 

(2tt)S 


The result is. 




r f * / w - mR m doj 

— I J U > r » ^ 

The velocity spectral density (per Hertz) is the Integrand: 
e^U 


Sd<“) . I . o.o) 


(110) 


If we set R « 0, we obtain the corresponding velocity 
spectral density in nonrotating coordinates: 

elU 


Sd(o) . j , 2. o.o) . 


(Ill) 


At this point, it is necessary to Introduce the integral 
scales of the turbulence. One such scale exists for each 
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of the three coordinate axes, but we need now consider 
only the axial and circumferential scales. They are de- 
fined as follows: 

I «(Xj, 0,0,0) dXj H Lj 

I ♦(O.XjjOjO) dXj 1 Lj . 

The velocity correlation function (j>(x,r,0) is plotted 
as a function of the axial coordinate xi in ;he upper curve 
of Fig, 7 • A plot of (j»(x,r,0) vs Xj would be qualitatively 
similar. Note that becomes very small when x, exceeds the 
axial integral scale Lj . The Fourier transform of is 
plotted in the lower part of Fig. 7 , where we see that the 
transform of the correlation function is small when kj 
exceeds 27 t,'‘L. . Similar behavior is demonstrated by the 
correlation function and its transform when plotted vs Xj 
and k,. In Eq. (Ill), the terms 

r.o) 

are small when either u>/U > 2 it/Li or m/r > 2 it/Lj . Con- 
sidered as a function of w, each term in the sum is 
similar in form, but of course becomes smaller as m/r 
increases. The result is a spectrum similar to the one 
shown in Fig. 6 . The spectrum decays monotonically as 
w Increases. The bandwidth of the spectrum is roughly 
“ 2itU/Li. When this is expressed as a fraction of the 
rotation rate, we obtain the following: 

'"'b . 2iT/n 

r rTTii 

The numerator is the period of rotation, whereas the de- 
nominator is roughly the time required for one coherent 
eddy to flow past the rotor. Thus, the bandwidth exceeds 
the rotation rate only if the axial length is small 
enough that the eddy passage time is less than the period 
of rotation. 

On the other hand, each term in the expression for 
the velocity spectrum in rotating coordinates [Eq. (110)] 
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peaks at a different multiple of the rotation rate. More- 
over, the mth term, considered as a function of frequency, 
is small whenever w is outside an interval of width 
4 ttU/L centered on mSl. But the frequency interval between 
peaks^is n, so the peaks are distinct whenever 4 ttU/Lj< 
or 


2TT/n ^ 1 

L^TU • 

Thus the spectrum seen in rotor-fixed coordinates consists 
of distinct peaks at multiples of the rotation frequency 
whenever the eddy passage time considerably exceeds the 
period of rotation. Equivalently, the spectrum in rotor- 
fixed coordinates is peaked whenever the bandwidth of the 
spectrum in nonrotating coordinates is appreciably less 
than the rotation rate. 

The heights of the successive peaks diminish as the 
rotation order m increases, and cease to be appreciable 
when m/r > 2ir/L2, or m > 2irr/L2. Many peaks occur only 
if the azimuthal Integral scale is considerably shorter 
than the rotor circumference, 2trr. The requirement for 
many distinct peaks is thus that the axial integral scale 
be long and the azimuthal integral scale be short. These 
conditions are met when atmospheric turbulence is drawn 
into the inlet of a turbofan operating at low or zero 
forward speed, as has been pointed out by Hanson (Ref. 

10) and others. The peaked spectrum which results is 
shown in Fig. 8 » when the axial Integral scale is 
short, a broad-band spectrum without sharp peaks is 
found . 

6.4 Spect»*al Density of Duct Modes Excited by Inlet 
Turbulence/Rotor Interaction 

In this section, we will derive an equation for the 
power spectrum of the sound generated by the rotor blades 
as they rotate through turbulence in the inlet duct. The 
starting point is Eq. (49), which gives the Fourier trans- 
form of the complex amplitude of the (ra,n)th duct acoustic 
mode generated by a given pressure distribution on the 
rotor blades. From this equation, we can write out an 
equation for the expected value of the (m,n)th mode (in 
the frequency domain): 
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<P„„(w) > ■ • 

n,m n,m 

rr_ fTr^ f+b* 

• j ^ R(r*,w) j ^ R(r”,v) I expCly(r* ,w)s' ] 


H 


H 


f+b" ( B-1 

• exp[-iv(r" ,u)) 2 "] •< I exp(i2TraiJ/B) 

lb" ^ 


B-1 

I exp(-i2TTinJl/B)<Ap. (r",2",u-mSJ) • 
i-0 ^ 

ApJ(r",s" ,v-ran) > I dz" dz' dr" dr’ , (112) 

where 

R(r,u)) i C (m/r)cosx-Yj^^j„(uj)sln xl (ll3) 

y(r,u)) = Y„ „(w) cosx+(m/r)slnx . (11^>) 

The expected value of the blade loading can be obtained 
from Eq. (55) ’• 

<Ap (r* ,z ' Ap*(r" ,z" ,v-mJl)> ■ 

J ^ 

(PflUj,)^ II f'(r’ , 2 * ,k,w-mn) exp(ljk«h) 
f*(r" , 2 " ,K,v-mJ5) exp(-lJlK‘h) • 

< w (r’,k,w-mn) w*(r" . (115) 

(2it)* (2Tt)» 
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The quantity in braces in Eq. (112). can be written out 
as follows: 

{Eq. (112)} = 

fill ,z’ ,k,(o-mfi) • 

f*(r” ,z” ,K,v-mfi) • 

B-1 

J exp[lj (2TTm/B+k*h) ] 

J=0 

B-1 

• J exp[iJ!,(27nn/B+K»h) ] • 


< w(r',k,u-mn) w*(r" ,K,v-mfl) > — ■ — . (116) 

(2tt)2 (27T)2 

The wavenumber-frequency spectrum [Eq. (106)] contains the 
factor 6(K-k) 6(kj-s/r), so in the expression above, we can 
set K = k^and kj * s/r. By noting that h * (2irr/B)g2> the 
summations over the indices J and i can be carried out 
explicitly; 


B-1 

\ expCi2ir(m+s) J/B] = B if m+s = pBj 
J*0 “ 0 otherwise 

where p is any integer. A similar result holds for the 
sun: over i. Thus, we can set s » pB - m in Eq. (106) and 
change the sum over all integers s to a sum over all 
Integers p. If we now Integrate over k and K, we obtain 
the following result: 

{Eq. (115)} = 


eJU(p,Uj,)=B = 
r 


6 (v-w) 


I f(r‘,z',k ,w-mn) 
P 


f»(r",z”,k ,v-mJ2) ♦(k„) 


(117) 
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where is defined as follows; 

V - I w-pBfl pB-m\ , . , 

% \-f-~ * V“ / • (118 ) 

It would appear from Eq. (117) that the elemental load- 
ing functions f(r,z,»**) and f *(r” ,z" , • • • ) must be computed 
anew for each value of the integer p. But, if we assume 
that the rotor is lightly loaded, so that the flow relative 
to the rotor blades is roughly aligned with the blade chords, 
then this is not the case. In the Integral equation from 
which the elemental loading function f(r,z,»»») is computed 
[see Eq. the wavenumber appears in two roles. First, 

the chordwise component of ijg appears in the expression for 
the chordwise distribution or the inflow velocity normal to 
the chord, which is the inhomogeneous or forcing term in 
the integral equation. The nominal Inflow velocity relative 
to the blades is Uy » U +( 2 re 2 . If we assume that the blade 
is lightly loaded^ then Ur will be roughly parallel to the 
blade chord, so that the" chordwise component of kp is 


Hr * , 


which is independent of p. The wavenumber also appears 
in the kernel function Itself (see Appendix b ) as the 
inter-blade phase angle a, defined as follows: 


a = 



h 


» 2TT(p-m/B) . 


But the kernel function is a periodic function of the 
Inter-blade phase angle, and the period is 2ir. Thus, the 
term 2irp can be deleted from a without af acting the' 
kernel function. It is apparent, then, th it the elemental 
loading function f(r,z,..,) does not depend on the integer 
p, so we can rewrite Eq. (117) as follows: 


(Eq. (115)} 


eJU(PoU^)^B^ 


6 (v-w) 


• f (r» ,z» ,kj^,w-mf2) f»(r” ,z" ,kj , \>-mn) 


• I I (kp,r,Ar) , 

where ko is obtained by setting p ■ 0 in Eq, (118). 


(119) 
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The expected value of the (m.n)th duct mode may now 
be written as 

£2u(p„U)2b2 6(v-u)) 

<p (ti)) p* (v)> * 

Pmn^ Pmn^ 


f^D rrj 3 r+D ■ 

J R(r’,w) J ^ R(r",w) I exp[ly(r* ,w)z’ ] 


H 


H 


n,m 

+b* 

-b’ 


pb" 

-b” 

• f(r”,z”,k dz" dz’ 


r+b" 

exp[-iu(r” ,(jj)z''] f (r’ ,z 'k^ jW-mfl) • 

\ / 


^ I (|)(k ,r,Ar) dr" dr' . 


r t. ^ '"p 

P ^ 


UL20 ) 


At this point, it is necessary to express the auto- 
correlation function (j>(x,r,Ar) as the product of three 
functions, each one of which depends only upon a single 
variable. That is,» 


(})(x,r,Ar) = * 


( 121 ) 


The three Integral scales Li, Lj, and L may be functions 
of the radius r. Procedures by which the functions <})i, 

4 > 2 , (}ij. and the length scales L;^, Lj, and Lr may be chosen 
to fit turbulence data collected from transducers placed 
in the inlet duct are described in Ref. ( 1 K The Fourier 
transform of 4» is then 


i(k,r,Ar) » (k^Li) ^jCk^L^) 

\ r/ 


( 122 ) 


‘ J 


«The representation given by Eq. (121) is physically meaning, 
lul only where Lj << 2Trr. 
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so 


^ ^ (u)-pBJ2)Lj 






(pB-m)L, 


\n) ■ 




(123) 


And the expected value of the (m,n)th mode Is 

pJU^B^ 6(v-w) 

* * 

^mn *^mn iir^k^ 

n,m 

f«-D fD f^'=' 

I R(r’,w) J R*(r”,uj) J f (r' ,z ' ,ko ,w-i 


mfl) 


H 


H 


-b' 


f+D” 

>xpCly (r* ,u))z' ] dz' j f (r" ,z” ,k ^ ,v-mfl) 

lb" 


exp[-iu(r",w)z"] dz" • ^ (^) 

(w-pBR)Lj] A f(pB-m)L< 


• *r(^) I K 


IT 




dr” dr’ . (12^) 


This rather cumbersome equation can be simplified if we 
assume that the radial Integral scale, Lj., is very much 
shorter than the radius of the duct, rp. If this is the 
case, both r’ and r” may be replaced by r = i(r’+r"), 
everywhere but in the argument of The integration 

variables can be changed from r' and r” to r and Ar » y, 
and the integration over y can be approximated as 


J <l>j,(y/Lj.)dr » Lp . (125) 

^00 
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When 


this is done, we obtain the following result; 




PjU®B26(v-w) 

n,m 


b^RMr.w) lC^(r,u))|* • (it) 




"H 

A. 

I 

P 


(o)— pBQ )L^ 

A 

'(pB-m)Lj 

— 0 J 

^ 2 

r 


dr 


where the chordwise integral Cjjj^(r,w) is defined as 


f+b 

C„«(r,o)) = f(r, 2 ,k ,a)-m«) • 

mn J ® 

-b 


dz 


expCiy(r,a))z] ^ 

The mean-square value of the (.Ti,n)th mode is given 


<P > 
*^mn 




rj/ v*.T dv d(i) 
exp[i(v-u))t ] j:n , 


which is 


<p* > 

*^mn 


I —2 I ® b*R*(r,uj) • 

i 8irrk2 J 


Sirrk^ _ _ 

n,m rjj 




IV 


(w-pBSI)Lj 


U 


(pB-m)L, 


dr 


d(t) 

T 5 i 7 • 


( 126 ) 


( 127 ) 

by 

( 128 ) 


( 129 ) 


55 



URjGiP^AL J 

OF POOR QUALITY 


The spectral density of the (m,n)th mode is the Integrand 
of the Integration with respect to (i)/2ir; 






n,m 


-r 


b^R^(r,d)) • 



(r,u)) 


r 


(U^/U)2 • 



(oj-pBS)Lj 

U 


Ak 



(pB-m)!^ 

r 


dr 


(130) 


Recalling the discussion of Sec. 6.3, (|>i is small unless 
the following inequality is satisfied: 



^ 2ir/nB 
L^D" 


Thus, if the axial Integral scale Li is long enough so that 
the eddy passage time, Li/U, is considerably longer than 
the blade passage interval, 2ir/f2B, the modal spectrum will 
contain distinct peaks at multiples^of the blade passage 
rate. Now looking at the function (j) 2 , we see that the 
peaks exist only at those multiples of the blade passage 
rate which satisfy the following inequality: 


lpB-m| 



Thus, if the azimuthal Integral scale is small compared to 
the me n circumference of the rotor, many peaks will 
^ appear in the spectrum. On the other hand, if L 2 is large, 
few peaks appear, and a mode selection process occurs by 
which modes with a number of azimuthal nodes equal to 
some multiple of the number of blades are most strongly 
excited. Finally, if Li is small, a broadband spectrum 
will result. 

Recall that the spectral density of the inlet tur- 
bulence sensed by a blade-mounted transducer tends to peak 
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at multiples of the rotor rotation rate, whereas acoustic 
mode excitation peaks at multiples of the blade passage 
rate. On a multi-b laded rotor these peaks occur at much 
higher frequencies. _ 


6.5 Sound Power Flux Generated by Inlet Turbulence/Rotor 
Interaction 

The sound power flux within the duct is obtained 
from Eq. (34): 


S(o)) 




I I 


mn 


(131) 


where Gnin(u) is the quantity in bracks given in Eq. (35). If 
we Introduce Snin(w) as the sound power flux per mode, 

i.e. , 


it(ri-rj) 

S (u) = ■■ rfi - G 

mn' ' p„U ran 


then the total sound power 
m and n; 


(w) 

flux is obtained by summing over 


S(w) = II 

m n 


(133) 


In the computer program, both S(w) and SiunCt*^) have been 
divided by PocSr^ to make them dimensionless. Using 
the dimensionless variables defined in Appendix A, the 
following expressions are obtained: 


S (w) 
mn' 


S(oj) ^ j j 

P.c^D P«°kD 


(IS*) 


Where 
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3**M B^w/n) 


!sn^ 

p.o!r^ /!^!t * „ Q y 

n,m/ 


o-o-’D ± M k. 




*|Cnm(x,w/n)12 e 




b^(x) 


Lj,(x) Lj(x) L,(x) • 

• I *|s - ") 


(pB-m)L, 


X 


dx 

X 


( 135 ) 


The chordwlse integral Cjjjn(x,u/n) is defined as 
follows : 


r+1 

C(x,w/n) » f(x,z/b,k,w) » 
-1 


exp[io^H(Yj^^j^cosx+| sinx)z/b] dz/b . ( 136 ) 

The elemental loading function f (...), the chordwlse Inte- 
gral Cjnn (...), and the modal power flux Smo/PjC^rd are all 
computed numerically. The procedures used are described 
in Chapter 8 . 
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CHAPTER 7 

ROTOR WAKE TURBULENCE 

7.1 Rotor Wake Turbulence Model 

Figure 9 shows a cylindrical slice of the rotor and 
stator opened out to form two linear cascades. Axes (xi,X 2 ) 
are fixed in the stator, while axes (Xi,X 2 ) move with the 
rotor blades. The equation relating the two sets of axes 
is 


X * X *K D + Orte, , (137) 

where C is a vector extending from the reference blade on 
the rotor to the reference vane on the stator. The re- 
ference blade and vane can be chosen arbitrarily, and D is 
defined only to within an arbitrary additive constant 
vector parallel to the X 2 and X 2 axes. In ^ig. 9 the 
centerlines of the rotor blades wakes are , cwn as dashed 
lines. These strike an angle x with the axis of rotation, 
as seen by an observer moving with the rotor blades. The 
unit vector o is normal to the wake centerlines, so the 
components of n along the Xi and X 2 axes are (-slnx, cosx). 
Note that X and n have- here different definitions than 
they had in Chapter 4, where they were defined with regard 
to the rotor blades themselves rather thar the rotor blade 
wakes. The two sets of definitions become identical only 
when the rotor is operating under ro load, so that the 
wakes are lined up with the blade chords. Since the equa- 
tions to be derived in this chapter are self-contained, 
no ambiguity should arise. 

The wakes are supposed to be convected with the 
nominal flow downstream of the rotor. Reynolds, et 
at. have conducted surveys of the turbulence within the 
wakes, using transducers fixed with respect to the rotor 
(Refs. 13,14) . These measurements indicate that the 
turbulence Intensity is maximum on the wake centerline, 
and diminishes monotonically with distance on either side 
of the wake centerline. The surveys were conducted 
relatively close to the Made trailing edges and show 
marked asymmetry of the wake widths on either side of the 
centerline. But at larger distances downstream, at 
distances comparable to the blade chord, the wakes become 
more syirimetric. 

Our mathematical model of rotor wake turbulence may 
be described as follows. Samples of the fluid velocity 
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FIGURE 9 . Rotor/SUtor Geometry for Wake Turbulence 
Analysis 
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downstream of the rotor, as seen in a rotating frame of 
reference, are assumed to have the following form: 

w(r,X,t) » WP(X«n) g(r,X-Wt) . ^ 

In this equation, w(r,X,t) is the component of the fluc- 
tuating velocity downstream of the rotor which is normal 
to the stator vanes, while W is the magnitude of the mean 
flow velocity, and eyr is the rms turbulence intensity 
on the wake centerline. The function g(r,X-Wt) is a 
random function of the variables indicated^ and P(X‘n) 
is a deterministic function chosen so as to model the 
variation of the turbulence intensity with distance 
from the wake centerlines. Note that X»n is constant on 
planes parallel to the wake centerlines.** To sustain 
the definition of Cw as the rms turbulence intensity, 
P(JC»r}) must be normalized so that P(0) » 1. Because 
all the wakes are assumed to be identical, P(X*D) must 
be a periodic function, the period being the dist inoe 
between wake centerlines, measured normal to the center- 
lines themselves. 

It is convenient to express F as an infinite series 
of identical ”haystack” functions, each representing 
the variation of the turbulence intensity across one of 
the wakes. The normal distance between wakes is 
2ur cosx/B, so P(X*q) can be written as follows: 


P(X*n) 

•V 



m 


( X*n-mh cosx \ 
6cosx j 


(139) 


where h = 2Trr/B. The denominator of the argument of 
f ( ), 6cosx, is the width of the wake turbulence intensity, 
measured normal to the wake centerlines. Thus, 6 itself 
is the corresponding width measured azlmuthally. One 
example of a candidate "haystack” function, and in fact 
the one currently being used in the computer program, is 
the well-known Gaussian profile, 

f(y/5) * exp[-ir(y/6)^] , ( 

where y « X«n. As defined here, 6 is equivalent to an 
Integral scale, in the sense that 
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f(y/6) dy • 5 . 


( 1 ^ 1 ) 


Because F(X»q) Is a periodic function, it can be 
expressed as a Fourier series, which can be reduced to 
the form shown here: 

■ hf I ‘(^) J-?) . (I''?) 

where f(a) is the Fourier transform of the "haystack" 
function f(x): 

f(a) i ff(x) exp(-iax) dx . (1^3) 


For example, if f(x) ■ exp(-irx*), then f(o) ■ exp(-o®/4Tr) . 

The random function g(r,X-Wt) is assumed to be a 
stationary random function, with respect to the spatial 
variable X* Therefore, the correlation function is a 
function only of X-Y, and not of X and Y separately. 

Thus, can be written as follows; 

<wCr^,X-Wt) w(rj,Y-Wx)> - ^^[x-Y-W(t-T ) , Ar] , 

where r ■ i(ri+r 2 ). Once again, the duct is annular, 
so must be a periodic function, Just as ♦jj was in 
Chapter 6 . ♦y, also can be written as an Fourier series, 

just as was. Thus, we have 


The Fourier coefficients ♦ are proportional to the Fourier 
transform of the nonperlodTc function 4> U,Ar), aeflned as 
follows, 


^(X,Ar) • ♦ (X,Ar) 

W 

if 

1 3Ci 1 <irr 

W) 

• 0 

if 

|Xil> r 


That is, 







(li*7) 
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' ■?ir ft »"'/J^»^r)®xp(lmX2/r) (1^<8) 

and the Fourier transform of t is 

♦^(A,Ar) = -iz$(X,Ar) 6(X2-m/r). (1^9) 


7.2 Wavenumber-Frequency Spectrum in Stator-Fixed Coordinates 


In order to calculate the power spectral density of 
the sound generated hy the interaction of rotor wake turbu- 
lence with the stator vanes, we must calculate the expected 
value of the wake wavenumber-frequency spectrum in a 
nonrotating reference frame: 


<w(r^ 


A 


,k,u)) w»(r,,K,v) 

^ * *v 


> = 



<w(r ,X-Wt) w(r,,Y-Wx) 

X ^ 2 -A* 


( 150 ) 


exp(la)t-lvT-lk*X+iK*Y) d^x d^y didt 


The expected value of the wake turbulence, as discussed in 
the previous section, is 


<w(r.,X-Wt) w*(r,,Y-Wx)> = 
e,y^F(X*..) P«(Y*n) • 


* • ^151) 


Prom here on, we will assume that the radial integral scale 
of the wake turbulence is small enough so that we can set 
ri ■ r 2 ■ r except in the last argument of above. To 

reduce the Integrals in £q. (150), substitute Eq.(151) in for 
the Integrand, and change the integration variables to X and 
Y, where 

X ■ X + D + f<rte. 


Y - y + D + nrxe, . (152) 

«•**«•* ^2 

Having done this, introduce a second change of variables, 
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5 « X - Y - W(t-T) , (153) 

^ ^ ^ 

where ^ replaces X. The integration with respect to t and t 
can now be carried out explicitly by recognizing the integral 
expression for the delta function [Eq. (98)]. The result at 
this point is 


<w(r-,k,a)) w(r.,K,v)> = 

(2tt) 2 5(o)-W*k+nrk,) 6(v-W»k+QrK, ) • 

• nil •n] F*(y*n) *(0 • 

exp[i(K-k)‘Y-ik»CJ dH d^Y • 

^ 


• exp[i(K-k)*D] 

*>• M 


( 154 ) 


The remaining inte^,- ations can be completed by substituting 
for F and F* the Fourier series as given in Eq. (142). The 
result, using Eq. (l49) for the Fourier transform of is 
shown below: 


<w(rj,k,u)) w*(r^,K,v)' 


6 (a)-W*k+Ork, ) 6 (v-W*k+nrK, ) 

«w ^ 2 <** Z 



(155) 
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Finally, it is worthwhile to introduce Up, the mean or 
nominal velocity relative to the stator^ 

U„ = W - Sire, (156) 

and to re-arrange the arguments of the five delta functions 
into a form more convenient for subsequent calculations. 

The final result is shown below: 


<w(r^,k,(u) w(rj,K,v)> 


(2ire„WB6)2 


w 


6(w-u *k) • 


exp|^~9^D*n 
'^•rcosx - - 


[h- Mr s] I «('' - 


(157) 


7.3 Wake Turbulence Velocity Spectral Density in 
Rotating and Nonrotating Coordinates 

In Sec. 6.3, the velocity spectral densities of inlet 
turbulence in fixed and rotating reference frames were 
derived from the inlet turbulence model adopted, in order 
to discuss the effects of rotation on the spectra. In 
this section, the velocity spectral densities of rotor wake 
turbulence also will be calculated in fixed and rotating 
coordinates, for purposes of comparison. 

The mean-square velocity of the wake turbulence can 
be calculated by setting ri ® r ‘2 ■ r in Eq. (157), multi- 
plying it by (2ir)* exp[l(v-aj)t+l(k-K) »X] , and integrating 
over V, K, and k: 
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(158) 




<w*(r,x,t)> = 

llllf <w(r,k,u) w*(r,K,v)> • 


exp[l(v-w)t+i(k-K)«X] 

' ^ " (2tt)® 


The result is 

<w^(r,x,t)> * 


(2iTr)®Uj,j^ m 




exp|i^S-^ (x+D)«n + in(il-m) Bt[ 
( rcosx f 


I 



mB 


rcosx I 


n\ dw 

n.Ol 


3 


( 159 ) 


where the wavenumber ks Is defined by the following two 
equations; 


— s 


U) 


e *k 
Z2 tTs 


s+mB 

r 


( 160 ) 


Behind the rotor, the mean-square velocity at a point fixed 
in the duct fluctuates at multiples of the blade passage 
rate, JJB. The time-averaged mean-square velocity is 
obtained by setting 1 « m and deleting the summation over 
Jt. If we delete the integration with respect to u, we 
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obtain the time-averaged rotor wake velocity spectral 
density. In the expression below, the components of kg 
have been written out explicitly. 


(2irr) 


I I 
s m 


^/mB6\P * / to+mBJ2 ^r2 s. s.r.O 


( 161 ) 


To simplify matters, assume that the .rotor is operating at 
no load, so that the azimuthal velocity in the wake, as 
seen by an observer fixed in the duct, is zero. Then the 
first argument of ^ depends only on the index m: 




The spectrjp consists of the sum of several "haystack" 
functions (j)i (w+mBn/Uj. i. . . ) centered at integer multiples of 
the blade passage rate. The width of these "haystacks" 
is roughly 4ir Un/Li, where Li is the axial Integral scale 
of the turbulence, and the frequency Increment between them 
is QB, so the peaks are distinct only if 



The first factor on the right-hand side of this in- 
equality is 27rr/B ■ h, the gap between the rotor blades. 


( 162 ) 


nRlGSNAl. PA3E l-a 


and the second is twice the ratio of the axial velocity be- 
hind the rotor to the rotational tip speed of the rotor. 

This latter ratio is roughly one half, so the peaks are dis- 
tinct only if the Integral scale of the wake turbulence is 
greater than the rotor blade gap. But ^he amplitude of 
the mth "haystack" is proportional to lf(mB6/r)l*, which 
is small unless mB6/r < 2tt, or m < 2iTr/B(S = h/d, where 
h * 2TTr'B is tb •) blade gap. This means that the high- 
order peaks are attenuated unless the width of the wakes, 

6, is much smaller than the blade gap. Therefore, if we 
assume that the integral scales of the turbulence within 
the wakes are comparable to the wake width, then the 
velocity spectral density, in a nonrotating frame of 
reference, will not contain pronounced peaks at multiples 
of the blade passage rate. On the other hand, if the 
wake width is considerably smaller than the rotor blade 
gap, but the axial length scale is larger, then a "peaky" 
spectrum will result. 

7.4 Spectral Density of Duct Modes Excited by 
Rotor Wake/Stator Interaction 

In this section, we will derive an equation for the 
modal spectral density of the duct modes excited by the 
interaction of the rotor blade wakes and the stator 
vanes. The starting point is Eq. (50), which gives the 
Fourier transform of the complex amplitude of the (m,n)th 
duct mode generated by a given pressure distribution on 
the stator vanes. Using this equation, we can write out 
an equation for the expected value (or ensemble average) 
of the (m,n)th mode in the frequency domain: 



(a,) p»^(v)> 



n,m n,m 


r+b ’ 

R(r",v) I exp[iy(r’ ,o))z’] 
-b' 


( V-1 

",o))z"] • < \ exp(i2TmJ/V) 

( J=0 


(Cont. ) 
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V-1 

* exp(i 2 TOJl/V) < Pj(r',2’,o3) • 


’ Pj(r«,2'',v) >1 dz” dz' dr" dr* 


(163) 


and y(r,w) are defined in Eqs. (116,117) * The 
( 54 “ '=^" be’obtiined 


<Apj(r', 2 *u) Apf (r'',z",v)> = 


(Po^r.^‘ il ^’(r',z',k,u) exp(ijk*H) 
• II f*(r",2",K,v) exp(-i)lK*H) • 


- £. 

• <w(r’,k,u3) w«(r",K,v)> 


d^k d^K 
(27T)'* 


(16A) 


and the quantity in braces in Ea 
out as follows: 


( 163 ) can be written 


{Eq. (163)} * 

^Po^r^' III! f(r’,z»,k, ) f«(r",z”,K, ) • 

V-1 

I exp[i,l (27rin/V+k *H)] 
j»0 ~ - 

V-1 

I expClZ(2irm/V+K»H)] • 

Jl»0 

<w(r’,k,w) w*(r",K,v) > 

(2U)** • (165) 

^But substitute - 0 for x. 

i 
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The complex wavenumber-frequency spectrum was derived 
in Sec. 5.2; it is rewritten below with different summatj n 
indices to avoid confusion later on. 


<w(r’,k,u)) w*(r”,K,v)> = 


(2ire„WB6)2 


w 


6(w-U„‘k) • 

•vI* *V 


_ [Km -1,) 


k - 


m^B \ 
rcosx 


6 C v-o)-n ( Jl -m ) B ] 


k-K- 


rcosx 


■h I 6 


/ s+m B\ 

(k- — ) • 


(166) 


Equation (165 contains two summations over the number 
of stator vanes, V. These can be summed explicitly if we 
note that the wavenumber-frequency spectrum contains the 
factor 


6[k-K-(m.-Jl, ) Bn/rcosx] 



s+mB 




which means that we can set ki - (s+mi3)/r and Kj * (s+)liB)/r 
everywhere they appear. Also, the intervene gap on the 
stator is H = (27rrg2/V), so 


k*H » 2iT(s+m B)/V 

^ X 


K»H » 27t(s+JI,B)/V 


and 
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V-1 

I exp[iJ(2Trm/V+k-H)] 

J=0 

= Vifm + s+ mB»pV 

1 *^i 

* 0 otherwise; and 


V-1 

J exp[li(2mn/V+K*H)] 

A=0 

* V If m + s + HjB = q^V 
» 0 otherwise, 

where px and qi are arbitrary integers, positive, negative, 
or zero. If we eliminate m and s from the equations 

m + s + m^^B = p^V 

m + s + )lj^B = q^V , 


we obtain the following relationship between m^, and 
Pi> ‘li* 


(m -i. )B = (p -q )V . 
11 11 


This equation has integer solutions for arbitary integers 
mi and Ai only if pi - m * JB for some integer J. Then 
mi - Ai = JV. If we set ii = mi - jV, then s = qiV - miB 
- JVB - m. The summations over l\ ana si can thus be 
changed into summations over qi and J. When this is done, 
the quantity in braces in Eq.(l63) becomes 


{Eq. (163)} 


(2Trr)^r 


l.l 

"»i ^1 


I 



■(m.-JV)B6' 

\~/ ' 

r 


(Cont, /' 
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I f (r ' ,z' ,k,w) II f*(r'',z”,K,v) • 

/ m B \ r (q V-JBV-m) 

‘ ? 

«(?-!; * 5 ^ 5 ) d^K 


6(u)-JBV-U ‘k) d^K d^k 

•>i»r 'w 


6(v-cj-JBV) 


• exp^i 


n.D 

rcosx - - 


Those terms in this equation with J “ 0 represent fluctua- 
tions at frequencies JBVn. The time-averaged modal 
spectral density is obtained by eliminating the summation 
on J and setting J » 0. Because of the factor 6(K-k) 
6(v-oj), we can set B = if and v * a) wherever they appear. 
Integrating with respect to K and k, we obtain 

{Eq. (163)}* — 6(v-o)) • 

(2iTr)^r 


A / m, 

I I fU 

mi qi \ 


m,B6\ 


• f (r* ,z* ,k,u)) 


m.B \ 

f»(r",z",k,u) (t»|k - n,Arj , 


(167) 
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where nov Js is a specific vector whose components in the 
axial and circumferential directions are 


* V V I 

^2 “ (qjV-m)/r . (l68) 


Equation (l67)-can be considerably simplified if we 
accept two approximations. The first is t,hat the radial 
Integral scale of the wake turbulence is much smaller 
than the stator vane chord. If this is so. we can 
replace r* and r” by their mean value r = i(r’+r'’) 
everywhere but in the last argument of the correlation 
function <|>. The second approximation is that the stator 
vanes are roughly aligned with the mean velocity in the 
wakes, so that the chordwise component of the wavenumber 
can be approximated as 


= oi/Uj. . 


This expression contains neither of the summation indices 
mi or Qi. The inter-blade phase angle, which is the 
product of the circumferential component of the wavenumber 
and the inter-blade gap, does contain Qj, but only as an 
Integral multiple of 2 tt. But the kernel function Kc is 
a periodic function of the inter-blade phase angle, and 
the period is 2ir, The Indlclal chordwise loading 
function f(r,z,k,u)) is therefore independent of m, and 
Qi, and so can be removed from the summations in Eq. 

(167). We then have 


{Eq. (163)} 


I Jir / 


6 (v-ui) 
-rU— 


(Cont. ) 
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m,B 


n,Ar) 


I 1 ~ rcosx ~ 

mi Qi 

expected value of into 

be written as shown below, by substituting ^q. ^ 

Eq. (163): 






. 2^ 

/ 4r^kJ^„(u>) 


■^D 

rr^ Uj(r)W(r) 

° 4'«(Xm * 

'U_ (r)r’ " ’ 

•H 


(| cos0+Yn^mS-^n0)^ • 

(r)l^ • 

I I 

m, q. 


dr' dr 


In this eauation, Cnm(x*) is the chordwise integral of the 
elemental ‘loading function, 


c (r) - f(r,z,k,w) exp(iyz) dz/b 

mn J 

-b 


(171) 


and 


u = 


Y COS0 - ? sln0 
’n,ra r 


(172) 
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Ir. Chap. 6, the correlatl-in function of the inlet 
turbulence was approximated as the product of three func- 
tions, each of which depended on only one variable. The 
object of this procedure is to facilitate the use of 
experimental data to select the appropriate correlation 
function. We will adopt the same tactic here: 


(|)(x,Ar) a ^j,(Ar/L^) (173) 

so that 

^ /s j\ 

4>(:.,Ar) a LjL 2 <|.i(kjL^) .^^(Ar/L^) (^a) 


and 


s m Bn 

( 17>lb ) 


where 


= U IT 

ri r 


2 


q^V-m^B-m 


(175) 


By changing the integration variables in Eq. (170) from r' 
and r" to r = J(r'+r") and Ar, and by approximating the 
Integral over Ar by the following infinite integral. 


4>(Ar/Lp) d(Ar) = 


the double Integration over the radius in Eq. (170) can be 
reduced to a single integral: 


<p (w) p* (v)> a 
^mn' , mr. ' 

(C-'nt, ) 


75 



I 


r 


ORIGINAL PAGE IS 
OF POOR QUALITY 




■ L L L 

_(0)) ^ ^ 
n,m 


U^(r)W(r) 

ijj u^^(r)r® *"»" 


i- oos6+Yn_„sme)' b=(r) |0^(r)|^ 
l^,m B5\|^ « 

^ ^ r "7“ ^ • 

' \ / ' 


In this equation, the chordwlse integral Cnm(r,( 
just as it was in Chapter 6: 


r+b 

Cmn(i'»w) * I f(r,z,k,ai) exp(iyz) dz/b , 

-b 


where now 


y * Y cos 6 - — sin0 
'n,m r 

The mean-squared value of the (ni,n)th mode 
calculated by integrating ov, • u and v: 

Fa a 

exp[l(v-(u)t] || . 

The result is 


( 176 ) 

i) is defined 

(177) 

( 178 ) 
can be 

(179) 
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U^(r)W(r) 

i ^(r)r* 


tp^(lC_ r>^) 
m,n 


|| COS0 + Yn^jjj(w)sinej^ (r) | Cj,j^(r,w) | 

K m,B6\I^ ^ 


d(i) 

2ir 


( 180 ) 


and the modal spectral density cf the (m,n)th mode is the 
integrand of the integration with respect to a)/27r: 


mn 


\ I 


I 


rn U^(r)W(r) 


(r)r^ 


H 


ri 


COS0 + b^(r) |Cjj^^(r,w) 1 ^ 

• I I KH— ♦jCX^L^) dr . USD 

m^ ' \ /' 


As discussed in the previous section, the modal spectral 
density will be broadband in nature unless the rotor blade 
wake width is less than the rotor blade spacing, and the 
axial Integral scale of the wake turbulence is considerably 
larger than the rotor blade gap. 
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7.5 Sound Power Flux Generated by the Interaction Between 
the Rotor Blade Wake Turbulence and the 
Stator Vanes 


The sound power flux within the duct is obtained from 
Eq. ( 3 ii): 


7T ( 3? ^ ^ ^ ^ 

b(u)) = ? I I Gnin^“^ » 

in n 




(182) 


where Gj^pCcj) is the quantity in braces in Eq. (35)* If* 
we introduce Snm(w) as the sound power flux per mode, 

l.e. , 


Tr(rl-rA) 

S (w) = -— 7 ? — <3 (w) P_„(u) 

mn p u mn mn 


(183) 


then the total sound power flux is obtained by summing 
over m and n: 


S(w) » I I 
m n 


(184) 


In the computer program, both S(aj) and Snm(‘*>) have been 
divided by poc^rj to make them dimensionless. Using 
the dimensionless variables defined in Appendix A , the 
following expressions are obtained: 


S(gj) , 

P c^r^ 

0 0 D 


+00 00 

= I I 

m*-« n=l 




(185) 


where 




(? ( o)/fl ) ( £^BV5 ) 

-^±Mlc I 

w n,m I 


( Cont^ ) 


(136) 
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dx 




The numerical methods used in compute this integral are 
discussed in Chapter 8 . 
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COMPUTER PROGRAMS OF POOR QUALITY 

Computer programs have been written which compute the 
i sound power per mode and the total sound power flux up- 

stream and downstream of a turbofan for the three noise 
K source mechanisms discussed in Chapters 5,6, and 7 . These 

computer programs are discussed in general terms in this 
' Chapter; FORTRAN listings and sample displays of the 

•; input and output data are provided in Volume 2. 

t 8.1 Rotor Mean Wake Velocity Oeficit 

The mean wake program computes the complex amplitude 
of each propagating mode excited by the interaction of the 
stator vanes with the mean velocity deficit wake of the 
' rotor. The sound power flux per mode is also computed, 

and by summing over alj. propagating modes, the total 
^ sound power flux is obtained. Recall that a multivaned 

stator excites only a subset of the propagating modes at 
any given frequency; specifically, only modes whose number 
of diametral nodes m is related to the number of rotor 
blades and stator vanes by the equation m = pV - qB, 
where p and q are arbitrary Integers. The program takes 
this selection mechanism into account in choosing which 
mode amplitudes to compute. 

A dimensionless version of Eq. (87) is used to 
compute the complex amplitude of each of the modes at the 
first three harmonics of the blade passage rate. The 
dimensionless form of Eq. (87) is written out below: 




(Uj./U)2 V 


27r(l-o2) k„ „ „ 
r n,m,q 


[ 'I' (X x)/— COJ3 6 + Y sin 9^ • 

J ^m^ mn ' ^x ’n,m,q 


exp 


i Y 6 

’n,m,q 1 x 2 




(187) 


Given the modal amplitudes Pmnq/PU*, the tota^. sound power 
is obtained from Eq. (28). In dlmer.cionless form, this is 


/> 
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Power _ j 
m 




m,n,q 


( 188 ) 


where Gnm(qBfi) is the quantity in braces in Eq. (28) or 
(29). 


>1 

I 


I 

a 

?< 


-i 

I 


In Eq. (187), Cnmq(x) is the chortiwise integral of 
the elemental blade loading function fq(r,z) multiplied by 
an exponential function; 


C 


mnq 


(x) 



fq(r,s) exp 


i 


Oc'o 

2 


'n,m,q 


cos8 - 7 sine 


) 


2 

F 


dz 


Because fq(r,z) has a square root singularity at the leading 
edge (z--b), it is necessary to change the integration 
variable from z to \(i, where z = b costii. We then have 




-r 


f,(r,z) 


sintl; 


* 


o.b 

/ \ 1 

* 

exp 

^ T" 

\^n,m,q ''p®® “ I sin0jcos.|;J 


di|; 


(189) 


This integral is easy to compute numerically using Simpson’s 
rule, because the product fq(r ,z)sinii) is finite at the 
leading edge (\J/»0,. The spanwise integration requires 
special treatment, however, because Wq is the qth Fourier 
coefficient of the component of the rotor mean wa!:e normal 
to the stator chord. 


Wq « Wq sin(0+x) exp - 

contains a phase angle which varies rapidly over the span 
when the separation between rotor and stator is large or 
when the rotor wake angle x is a strong function of the 
radius r. To handle this rapidly varying phase, an 


; 

TV 

'd 



u 
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adaptation of Filon's rule of integration has been used 
in place of Simpson’s rule for the spanwise integration 
(see Appendix G). 

8.2 Inlet Turbulence and Rotor Wake Turbulence 

Random loads on the rotor or stator blades caused 
by convected turbulence generate continuous noise spectra, 
as opposed to the purely tonal noise generated by the 
rotor mean velocity deficit wakes. To plot these con- 
tinuous spectra accurately, it is necessary to compute 
the modal power flux In the turbofan duct at many more 
frequencies than are required for the mean wake tone 
noise. To save computer time, therefore, it was decided 
first vO compute and store the elemental blade loading 
functions f(r,z,k,o)) over a pr-'determined range of 
frequency and wavenumber, and then to interpolate between 
these values to compute the power spectral density at 
intermediate values of the frequency. The advantage 
gained by doing this arises from the fact that the elemental 
loading functions can be computed and stored once and for 
all for a given rotor /stator geometry. Once this has been 
accomplished, parameter variations (frequency, turbulence 
length scales, etc.) for the chosen rotor/stator geometry 
can be conducted conveniently and economically. The 
elemental loading functions have been found to be reason- 
ably smoothly varying functions of the fi’equency (and 
for convected turbulence frequency and wavenumber are 
related), so acceptable accuracy (within 3 %) can be 
obtained by computing and storing the loading functions 
at integer multiples of the rotor rotation rate fi. Care 
must be taken, however, to assure that the array of 
stored loading functions encompasses the required range of 
frequency, especially for the inlet turbulence case, 
because in that situation the frequency appears in the 
combination <o - mfl, where m, the number of diametral 
nodes, can take on large positive and negative values. 

The inlet turbulence and rotor wake turbulence pro- 
grams print out the spectral density of the power flux 
per mode for each propagating mode, and by summation, the 
spectral density of the total sound power flux. Dimen- 
sionless expressions I'or the spectral density of the power 
flux per mode are given in Eq. (135) (inlet turbulence) 
and Eq. (l86) (rotor wake turbulence). Both equations 
require numerical Integration over the chord and span 
of the blade (or vane). The integral over the chord is 
basically the same as the integral required in the mean 
wake program; the same change of integral variable is used 
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to eliminate the singularity at the leading edge. Once 
this change of variable has been effected, a straightforward 
application of Simpson's rule is sufficient for both the 
chordwise and spanwlse Integrations. 

8.3 Integral Equation for the Blade Loading 

As discussed in Chapter 3, the elemental blade loading 
function f ( • ,y/b, . . . )* is the solution of an inteLral 
equation, 

f +b 

exp(ik^z) = J K^(«, f(*,y/b,**«) dy/b , 

-b 


whose kernel function Kc( * , (z-y ) (b, • • • ) is derived in 
Appendix B. Given the elemental blade loading function, 
the pressure on the blades is obtained by multiplying 

f(*,y/b,***) by poUr w(r,k,u)) and inverting the Fourier 
transforms with respect to time and/or space as appropriate. 

The integral equation contains two basic difficulties. 
They are 


a. the solution f(»,z/b,**0 is singular as (z+b)”^ 
at the leading edge U=-b), and 

b. the kernel function Kr ( • ,y/b, • • • ) contains both 
a Cauchy singularity (l/y) and a logarithm:* c 
singularity at y = 0. 

Difficulty (a) is circumvented by introducing the 
independent variable transformations 


z/b = cos e 

y/b = cos 

« 

Equation (B.l) then becomes 


*Ir previous chapters, the elemental loading function f( ) 

has been written as having dimensional arguments. Use of 
dimensionless arguments in Chapter 8 is convenient and should 

:.-use no confusion. 
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expdk^b COS0) ® ( • ,cos9-cosi(<, • • • ) F(\|<) dt|) , 

“ (190) 

where F(i|^) ® f ( • ,cosit», • • • )sin4». Whereas f ( • ,cosip, • • • ) is 
singular at = it, Fd) is not. Thus the integral equation 
is solved for F(\j;) rather than for f(**0 itself. 

To solve Eq. (19C), the method of collocation is used. 
(Refs. 18, 19). That is, the integral on the right-hand 
side of Eq. ( 130) is required to equal the forcing function 
on the left-hand side at the N points. 

®m “ (m-l/i^jir ; m - 1,**»N , (191) 


The equations to be solved are then 


expdk^b cosSj^) = ~ 1,2,***N 


( 192 ) 


where 


1 ( 0 ) 


7T 

K (•,cos0-cosip,'**)F(ii;) diij 


0 

In computing 1(9), the left-hand side of Eq. 
points 


(193) 

( 19 3) , the 


'^n ~ i “ 0,1,2, •••N 


are used.* Having the collocation and integration points 
evenly interlaced has the dual effects of avoiding the 
point 4; = 0, where Kc is singular, and of Insuring that the 
solution obtained satisfies the Kutta condition, which 
requires that f(***) » 0 at the trailing edge. 


•Because F(0) ■ -, use of N -*• 1 integration points as 
called for in Eq. (193) Introduces only N unknowns f(fir/n), 
n’sl , . . .N . 
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The logarithmic portion of the kernel function must 
be singled out for special treatment. Thus, 


+ G(*,y/b,‘»*) log |y/b| 
and 1(6) * 1^(0) + 1^(6) , 

where 


ij(e) 


•Tr 

Kj ( • ,cos9-cosij;, • • • )F(t|)) d»|/ 

0 


1^(9) 


TT 

G( • ,cos9-cosi|^, • • • ) • 


log 1 COS0-COSij^ lF(l|)) dij; 

The trapezoidal rule suffices to compute 1^(0), 


N 

I (9) = I B K ( • ,cos6-cosil< , • • • )F(4^ ) u/N 

i ii 1 ii n 


where Bj, » Bj. = 1/2, while B = 1 otherwise; but a special 
integration rule, devised b^ Whitehead (Refs. l8, 19) 
is required for the logarithmic portion of the kernel. 
Whitehead's integration rule is given below, 

N 

1,(0) « I BpG( •,cos9-cosij<., • • • ) • 

^ 1=0 * 


log |cDs6-cosij^^I Sj^(9) tt/N , 

where the Bj^'s are the same weighting functions as defined 
above, and 
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Sj^(e) = - {iog(2) + 

N 

2 y (B /s) cos(s9) cos(s£tt/N)} 
s=l ® 

By combining 1^(0) and Il( 0), the integral equation is 
reduced to the ^following set of algebraic equations for 
the unknowns F(n-n/N): 


exp[l(k b) cos (m-1/2 )tt/N] = 

! F(n,./N) , 

n=o 

where 

\n * Bn^^c(j'*)+G(x)[S^(ejj^)-loglx| ] ir/N 

and 

X = cos(bjjj) - cosCij^j^) 


Any standard matrix package that handles equations with 
complex coefficients can be used to solve these equations. 

FORTRAN listings for computer suborograms which 
compute the elemental loading function f( ), based upon 
the equations set forth in this section are provided in 
Volume 2. These subroutines have been used to compute 
the blade lift and pioching moment coefficients for a 
cascade of oscillating blades in uniform flow, and for 
a rigid cascade of blades in a sinusoidal gust, for 
comparison with results of similar calculations published 
in the lite'^ature (Refs. 18, 19). Agreement was obtained 
to four sign, f leant figures in every case tried. Com- 
parisons have also been made with the loading function 
itself, published in graphical form by Ka.Ji (Ref. 5). 

Again agreement was excellent, within the limitations 
Imposed by the small scale in the figures provided in 
Kaji’s paper. 
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IriAPTER 9 

SAMPLE COMPUTATIONS 


No computations have been made with the Inlet Turbulence/ 

Rotor Interaction program, and only a few have been made with 
the Mean Wake/Stator Interaction and Wake Turbulence/Stator 
Interaction Programs. Typical results from these latter compu- 
tations are presented In this chapter only as examples of 
output data generated by the programs; a systematic set of 
parameter variation studies has not been performed. 

The Mean Wake/Stator Interaction Program computes the 
sound power per mode and the total sound power radiating up- 
stream and downstream of a turbofan at the blade passage 
frequency and Its first two harmonics. The sound radiating 
upstream from a 15 blade fan and 11 blade sta.or is displayed 
in Pig. lO’^. As indicated in the Figure, wh -eas at the blade 
passage frequency only one propagating mode is excited, many 
are excited at higher frequencies. 

Typical results obtained with the Wake Turbulence/Stator 
Interaction program are shown in Fig. 11. In this figure, 
the upstream sound power spectral density for various values 
of the turbulence length scales is plotted as a function of the 
frequency. The axial, radial and azimuthal length scales were 
assumed to be equal. The width of the turbulent wake is 
0 . 072 rr, and the rms turbulence intensity is assumed t' be 
12 of the nominal flow velocity. (The sound power is directly 
proportional to the rms turbulv-.nce intensity, so the sound power 
at any other intensity may be found by shifting the vertical 
scale in the figure.) It may be seen that increasing the turb- 
ulence length scales caused bands of noise to coalesce around 
harmonics of the blade passage frequncy. However, the length 
scales required to produce this effect are probably unrealistically 
long. Individual modal amplitudes are not plotted because of 
the large number of propagating modes (about 130 modes at three 
times blade passage frequency). 


•The fan/stator geometry is that of NASA Pan Stage 55 (Ref. 20). 
The input variables used for the computation may be found in 
Volume 2 (Mean Wake Program). 
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CONCLUDING REMARKS 


Equations have been derived for the amplitudes (or 
amplitude spectral densities) of the propagating duct modes 
excited by a turbofan (fan/stator stage) operating at 
subsonic tip speed within an infinite hard-walled annular 
duct. Three noise source mechanisms were considered: 

a) random noise generated by inlet turbulence 
impinging on the rotor blades; 

b) random noise generated by rotor blade wake 
turbulence impinging on the stator vanes; and 

c) tone noise generated by the mean velocity defect 
wakes of the rotor blades impinging on the stator 
vanes. 

Equations have also been derived for the acoustic pres- 
sure at an arbitrary point within the duct, and for the sound 
power flux (or sound power spectral density) in the duct 
either upstream of the fan or downstream of the stator. 

These latter equations require a summation over all those 
duct modes that are excited by the particular noise source 
mechanism of interest and which propagate at the selected 
freqvcacy . 

A package of FORTRAN computer programs has been 
developed which computes the duct mode amplitudes (or 
spectral densities of the mode amplitudes) for each of the 
three noise source mechanisms, and sums these mode 
amplitudes to compute the sound power flux (or sound power 
spectral density) within the duct. The inlet turbulence and 
rotor wake turbulence programs print out, at any desired 
frequency, the sound power spectral density per mode for all 
propagating modes and the total sound power spectral 
density. The mean wake program prints out, at the blade 
passage frequency and its first two harmonics, the mode 
amplitude and sound power per mode for each mode excited. 

The total sound power flux is also printed out at each 
frequency. 
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A very limited set of preliminary computations has been 
carried out using the mean wake and wake turbulence 
programs. The output data display the expected qualitative 
behavior, but no attempt has yet been made to make any 
systematic comparisons with experimental data. No 
computations have been made with the inlet turbulence 
program other than those required to verify that the program 
runs. 
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APPENDIX A 
LIST OF SYMBOLS 

mode amplitude (Eq. 23) 

# or rotor blades 

blade/vane semi-chord 

blade/vane semi-chord at radius r 


blade/vane semi-chord at tip 
chordwise Integral (Eq. 130) 
speed of sound 

displacement vector extending 
from rotor blade to stator 
vane (Eq. 62) 

axial distance from rotor to 
stator 

exponential function 

unit vectors (Fig. A) 

amplitude of turbulence In- 
tensity in rotor blade wake 

elemental blade/vane chord- 
wise pressure distribution 
(Eq. 55) 

elemental vane chordwise 
pressure distribution at qth 
harmonic of the blade passage 
frequency (Eq. 85 ) 

mean rotor wake velocity de- 
ficit profile (Chapter 3): 
or rotor wake rms turbulence 
intensity profile (Chapter 5) 


h 

4. 

' 'ii 

: 

i 

I 
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\ 


g (r, x-wt) 
H 


h 

I 

1 


See Eqs. (13^) and (35) 

See Eq. (l4l) 
inter-vane gap (stator) 
inter-blade gap (rotor) 
axial intensity 

/rr- 

mth order Bessel function of 
1st kind 


(uj) 

,m ^ 
k = (k^jk^) 



wavenumber vector 

See Eq. (39) 

wavenumber vector 

particular value of wavenumber 
vector as defined in Eq. (118> 


Turbulence length scale in the 
1th direction. 


M 


axial flow Mach No. 




Mach number of flow relative 
to blades/vanes 


Mn 


relative Mach number at the rotor 
blade tip. 


N 


unit vector normal to stator 
vanes 


n 


p ( w) 
^mn 


Pmn ^ ) 

tun 


R (r,w) 


unit vector normal to rotor 
blades (Chapter 4)j or normal 
to rotor wake centerlines 
(Chapter 3 and 5) 

acoustic pressure 

duct mode amplitude 

amplitude spectral density 
of (m,n)th mode 

See Eq. (116) 
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S(w) 
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- (Xi,X2) 
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radial coordinate 
hub radius 
duct radius 

total intensity spectral 
density 

intensity spectral density 
of (m,n)th mode 

time 

axial fluid velocity 

fluid velocity relative to 
rotor or stator vanes 

number of stator vanes 

velocity downstream of the 
rotor (defined in rotor-fixed 
coordinates) 

component of fluctuating 
velocity normal to rotor/stator 
blades 

rotor-fixed coordinates 

axial coordinate 

stator-fixed coordinates 

integration point (r,yj, $ ) 

mth order Bessel function of 
2nd kind 

chordwlse variable 


PAGS 

quality 


A- 3 
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z 
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GREEK LETTERS 


a 


6 

6 


r 


Y ( U) ) 

APs(r,y, . . . ) 


r 

6 


i » (61,62) 


e 


w 


e 

1 * (Xl,X2) 
y(r,oj) 

V 

i 

ir 


Fourier transform variable 
/l-M^ 

/l-Mr 

See Eq. 11 

pressure loading on blades 
(Eq. 53) 

TT irl . rg) 

half-velocity width of rotor 
mean wake (Chapter 5), or 
width of turbulent rotor wake 
(Chapter 7) 

blade/vane sweep 

rms Intensity of th>. component 
of the inlet turbulence normal 
to the blade surface. 

rms intensity of turbulence 
on cen..erlines of rotor wakes 
(component normal to stator 
vanes) 

stagger angle of stator vanes 

eigenvalues of annular duct 
modes 

Fourier transform variable 
See Eq. (Il7) 
radian frequency 
integration variable (Eq. 153) 

3.1^159. . . 
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nominal fluid density 
Interblade phase angle 

2b,j/rp 

time 

polar angle 

rotor blade stagger angle 

normal mode 

radian frequency 

rotor rotation rate 
(radlans/sec . ) 
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APPENDIX B 

KERNEL FUNCTION FOR A LINEAR CASCADE IN SUBSONIC FLOW 

The kernel function for a cascade of thin airfoils in 
oscillating subsonic flow has been dervied by several in- 
vestigators, using a variety of methods (Refs. 1-6). Of 
the various approaches employed to date, the method of 
Fourier transforms is perhaps the most straightforward. 

This procedure was used in Ref. 3, for example, but the 
inversion of the Fourier transform was accomplished numer- 
ically. More recently, Goldstein (Ref. 7) pointed out 
that the Fourier transform of the kernel function contains 
no branch points, so the inversion can be accomplished 
quite easily by using the Cauchy residue theorem. Gold- 
stein outlined the procedure to be followed, but did not 
actually carry out the calculation of the kernel function. 
The purpose of this Appendix is to set forth the details of 
the inversion, and to record the end result. For the con- 
venience of the reader, as well as to document the notation 
used, a bi‘ief derivation of the transform of the icernei 
function is also presented. 

The cascade geometry is shown in Pig. B.l. The air- 
foil semichord is b, and the gap between neighboring 
airfoils is bh, with components bhi projected along the 
chord, and bha normal to the chord. The airfoils are shown 
as having no camber, because the ultimate objective is 
to calculate the pressure field scattered by the cascade 
when it is subjected to vorticlty convected with the mean 
flow. However, other situations, such as a cascade of 
oscillating airfoils, can be handled as well. For the 
purpose of calculating the kernel function, we need only 
suppose that a known chordwise pressure distribution exists 
on each airfoil, and calculate the resulting velocity field. 

The first step is to calculate the upwash generated 
by a single airfoil. Let the pressure be given by the real 
part of p expdwt), and the corresponding velocity field 
be the real part of (wi,W 2 ) exp(iwt). Then p satisfies 
the convected wave equation, and the velocity field is 
related to the gradient of p through the momentum equations. 
These equations are* 


*A list of symbols for Appendix 3 3s provided on page 312. 


3-1 
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PIG. B.l- CASCADE GEOMETRY. 
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0 , 


(B.l) 


do. + Orsi;’ "J * ^ ° ■ '*' 

The following definition of the Fourier transform pair is 
chosen: 


P 



expC-iayj) 


dy, 


P 



exp(+iayj ) 


da 


(B.3) 


The Fourier transforms of Eq. (B.l)* and the second of 
Eqs. (B.2) are 


^ - X*p 0 , (B.4) 

dyj 

(ii»+loUj.) ^ ^ ' 0 . <B-5) 

where 

X s {a*-M^ ( a+w/Uj, ) * ) ^ . (B.6) 


The pressure p generated by an isolated airfoil is bounded 
at infinity, and, because the airfoils are assumed to have 
zero thickness, is antisymmetric in y^. A solution of tq. 
(B.4) satisfying these conditions is 


•All equation references are to equations in this Appendix. 
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p = _ A£Xoj_ expC-Xl^jl) • sgnCy^) , 


(B.7) 


provided that X is defined so that its real part is non- 
negative on the path of integration_used to invert the 
Fourier transform. In Eq. (B.7)j Ap(a' is the Fourier 
transform of the chordwise pressure distribution on the 
airfoil. (Ap is positive if the pressure is greatest on 
the lower face of the airfoil.) By eliminating p between 
Eqs. (B.5) and (B.7), we obtain the Fourier transform of 
the upwash generated by a single airfoil. The upwash 
itself is 


Wj(yj,y2) 


u. 



X Ap(g) 

21(o+u)/Uj,) PjUj 


• exp(iayj-X [yj 1 ) da 


The upwash generated by a cascade of airfoils, located at 
the points (yi,y,) = mCbhjjbhj), m = 0, ±1, ±2,..., is 
obtained by summing the contributions of the individual air 
foils: 


“r » I, ai(o+w/Uj.) 

• exp{ia(yj-mbhj ) - Xly^-mbh^j} da 

This infinite series can be summed analytically if 
the transformed pressure distributions on successive air- 
foils in the cascade [APm(V)] are related by a constant 
increment in phase angle. That is, for any integer m. 


Ap.(o) - Ap-(a) exp(ima) 
la ® 


9 
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where Apj(a) is the transform of the pressure on the "zeroth" 
or reference airfoil, and a (called the Interblade phase 
angle) is a constant, 0 _<o< 2 tt. To calculate the upwash near 
any selected blade, say blade s. Introduce 



and 


n = m - s 


Then 


^ (x,+sbhi,X2+sbh2) = 

expUso) f*" ‘ exp'laxi). 

" 27T J 0 2i(aWU_) 


-00 ^0 r 


(B.8) 


where the infinite series 


+eo 


S(a) = I exp{ino-A [x2-nbh2 1 -ianbhj} 


n=-® 


can be summed: 


S(a) - J 


^ j exp(|-A^-Ax2) . exp(jA_+Xx2) 

sinh(|A^) slnh(|A_) 


where 


A i ±Xh2 + l(o-abhj) . 


(B.9) 


(B.IO) 


Now, use the convolution theorem to calculate the Inverse of 
Eq. (B.8): 
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w. 


g- (Xj+sbhj jXj+sbhj) “ 

r 


f+b jjj 

exp(isa) K^(Xj-C,X 2 ) — 

ib 


Apo(C) 




(B, 


is the desired cascade kernel function: 
c 




^ /•+• XS exp ( lax j) 


-.00 


2i(a+w/Uj.) 


da 


(B, 


The factor b Is Included to make dimensionless. 

It is convenient to introduce dimensionless variables 
follows : 


wb /$2 Cj 

X 

III 

< X 



ab-M^K 

Xj i Xj/b 

• 

(B. 


The shifted transform variable a completes the square in A 
X - ^ ^a*-K» 

Then 


K 


c 


exp(lKMj,^j) 


+• yS exp(loXj) 
^ 2 (a+K/M^) 


da 


where now 


(B. 


Y = 


(o*-K*)^ 


(B. 


11 ) 


12 ) 


, as 


13 ) 


li*) 


15 ) 


B -6 
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[ exp(|-A^-6j.TX2) 
sinh(^A^) 

exp(jA_+Bj.YX2) 

sinh(^A_) 



(B.16) 


V 


Notice that yS is an even function of y, so that even 
though the integrand contains the variable 


y « , 

A 

there can be no branch points at a = ±K. Thus, if we apply 
the residue theorem to evaluate the kernel function by 
closing the path of integration on a large arc in the upper 
or lower 6i plane, no residual integrals around branch cuts 
appear; the kernel function is simply the sum of the resi- 
dues in the upper or lower half plane. To insure that the 
integrals on the arcs vanish as their radii are allowed to 
become infinitely large, the Integrand in Eq. (B.lM) must 
be modified. First, note that 

S « 1 ils 

so can be written as follows 


exp(iKMj,Xj) 

2 ,+• S exp(iaxj) ^ (B.17) 

. .2 ;; do 

3x* 2y(o+K/Mp) 


The Integral above can be evaluated via the residue 
theorem by using the contours shown in Pig. B.2, and the 
differentiation with respect to Xj carried out afterwards. 




I 

i ' 
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If Xi>0, the integral around the semicircle in the upper 
half plane vanishes when the radius of the contour goes 
to infinity. Thus, 


*0 ' - r 

r 

. /sum of residues in) 

3,2 \ upper half plane f * (B.18) 

2 

On the other hand, when Xi<0 the Integral around the contour 
in the lower half plane vanishes, so 


Kg = + ^ expUKM x.) . 

r 

3^ fsum of residues \ 

\in lower half plane/ (B.19) 

Tie integrand in Eq. (B.17) has poles at a = -K/M , 
and at points a^, where sinh(l/2A^) = 0: 

® 2mri (n * any integer). 


[c.f. Eq. (B.l6)]. When A+ is real, some of these poles lie 
directly on the real axis. To arrive at the correct expres- 
sion of Kg, it is necessary to invoke the casuality condition 
by stipulating that w has a small negative imaginary part. 
Once the kernel function has been evaluated, we can let 
Im(oj)-**0. The effect of this procedure is to eliminate the 
possibility that acoustic waves not generated by the cascade 
itself are inadvertantly Included in the solution. 

If lm(w)<0, the pole at a » -K/Mj, is clearly in the 
upper half plane. The residue is 


slnh(ejKh/M^)exp(-lKStj/M^) 


‘K 


^^r^ cosh(0*Khj/Mj,)-cos(r+Khj/Mp) 


(B.20) 


a 




I' 
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by 


The roots of the equation sinh(l/2A+) = 0 are given 


+ r h, 3^h, 
o- - -2-L 1 J 
n j» d 


(a=) 


Tji2lh 


(B.21) 


\ •• 


where 


i r - 2nir 
n 


d - /hJ+Bjh* 


(B.22) 


If the square root in Eq. (B.21) is defined as follows, 

[Re^®]^ 5 /R , 0 < 0 <2ir (B.23) 

where is the positive square root of R, then solutions 
having the plus (minus) sign in Eq. (B.21) are located in 
the upper (lower) half plane. In either case, the residue 
is 




2d" 


expda^Jj ) 




(B.2i|) 


Using Eqs. (B.l8) and (B.19), and carrying out the indi 
cated differentiations, we obtain the final result: 


B-10 
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2 v slnh 


*c * 


2M, 




2d' n-- 


/ WJKX.; 

l~) -x~) 

/B*Kh\ / Kh, \ 

.e„ 3 (r . 

)exp[l(a%KM^)xj ] 

(■: ■ W- • -i) 


If Xj > 0, and 


.2 


ejhj (a“ -K^)exp[i(o"+KM^)xj j 

K_ = — I 


2d* n*-® 




(B.25) 


(B.26) 


If Xj <0. Having calculated the kernal function for 
i;^(K) < 0, we can now let IM(K) ->■ 0. Then the roots 
become, using the branch of the square root given in a 
Eq. (B.22), " 



(B.27) 


if \T^/(i\ > K, and 



if |rn/dl < K. In Eqs. (B.27) and (B.28), is to be 
calculated using the upper__set of signs, end a^ using the 
lower set of signs. The / sign means the positive square 
root. 
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b 

d 

h. 


i 

K 

M 

r 

P 

S 

t 

(x, *Xj) 
(y^pyj) 


airfoil semi-chord 
speed of sound 
/hf+Bih? 

(inter-blade stagger gap)/b 
(inter-blade normal gap)/b 
/TI 

reduced frequency [Eq. (13)] 

Mach number 

pressure 

Eqs. (9, 16) 

time 

nominal fluid velocity 
perturbation fluid velocity 
coordinates [Eq. (13)] 
coordinates (Fig. 1) 



Fourier transform variable 

ob-M K 
r 

✓IImJ 


Eq. (10) 
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X 

7T 

Po 

a 

CO 


Eq. (6) 

3 . 1^4159 ... 
nominal flow density 
Interblade phase angle 
frequency 
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APPENDIX C 

NUMERICAL COMPUTATION OF THE NORMAL MODES IN AN 

ANNULAR DUCT 

As discussed in Chapter 2, the normal modes in an 
annular duct are functions of the form 

♦m' ® ’'n'W’ > ' = 


where Jm( ) and Y,„( ) are Bessel functions of the first 
and second kinds, and Kjnn are the roots of the following 
transcendental equation: 


i 

i 


•I 

4 


Y'(K„r„) 
m mn H m mn H 


m mn D m mn D 


= 0 


(C.2) 


The roots of Eq. (C.l) are found by first estimating 

the root as follows: 


Wd 




m 


I "^mjn- 


if n = 1 


^r^j + IT if n > 1 


(C.3) 


This estimate is refined by incrementing the estimated 
value of KmnrD by it/ 10 until the determlnent in Eq. (C.2) 
changes sign. The step size is then halved and changed 
in sign. This process continues until the absolute value 
of the determlnent is reduced to a preassigned value. 


Once the eigenvalue icmn has been computed, the con- 
stants A and B are assigned one of the following two sets 
of values: 


A 


B 


B. 1 


• i 



\ 


C-1 


A 



V- 
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Of these two sets of values, the one for which (A^+B^) is 
the smaller value is chosen. The desired normalization, 
namely , 

*’n< ' i 


is obtained by computing the value of the integral on the 
left-hand side of the equation above, using the formula 



The constants A and B are then divided by whatever common 
factor is required to produce the correct normalization. 


i 


I 


I 

i 


C-2 


i 
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APPUiDIX D 

FILON'S INTSGP.^TION RULE 

Filon’s integration rul* 'Applies to integrals of the 

^rm 


J f(x) sxp(igx) di (D.l) 

a 


wherein ^h^ phase ang‘. ... is large, so that use of the 
trapezoidal rule or i ^p .•.a's rule would require that the 
interval of inte gratis . ae divided into many subintervals 
to obtain an accurate answer. Filon's rule is obtained 
by assuming that the function f(x) can be approximated 
by a quadratic function, but the exponential function 
exp(igx) is integrated exactly, without approximation. 

To extend this procedure to Integrals wherein the phase 
angle is not a linear function of x, we need only approxi- 
mate the phase angle as a linear function of x within 
each subinterval. The most straightforward procedure 
is to approximate both f(x) and g(x) as linear functions 
of X. For example, to compute the following integral 

I = I f(x) expCir'c)] dx , 
a 


divide the Interva?. (a,b) into N equal subintervals of 
length h = (b-a)/N. The integral over the nth subinterval 
is then 

rx th 

» f(x) exp[ig(x)] dx . 


To compute this integral, approximate both f(x) and g(x) 
as linear functions: 


f(x) “ f^ + 

(‘•n*rfn> 

/»-v 

\ h / 

> 

g(x) » g„ + 

^®n+l"*n^ 


1 . 


D-1 
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where and so on. Carrying out the integrations, 

we obtain 


I„ . er.pdg,) , 


where 

So = 2 ^^n'*'®n+l^ 


^ ^ 2 ^®n+l"®n^ 

a = sinX ^ J sinX _ c osX \ 

\ 2 \^ 2X / 

and a* is the complex conjugate of a. When X 0, the 
integral over the nth subinterval becomes 

In “ I • 


which is the trapezoidal rule. 

If an even number of subintervals are used, it is also 
possible to approximate f(x) as a quadratic function (over 
any two neighboring subintervals) while leaving the phase 
angle g(x) as a linear function of x. Th*s Integration 
scheme reduces to Simpson's rule when g(x) is constant. 

Both the trapezoidal and the quadratic versions of Pillon’s 
rule have been tried; the quadratic version did not seem 
to Improve the convergence significantly over the simpler 
trapezoidal rule, so the latter is used in the rotor mean 
wake program. 


